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INTRODUCTION 

This report is a summary of the work performed under NASA grant NAS8-39716 (Liquid 
Metal Diffusion). This research was to advance the understanding of diffusion mechanisms in 
liquid metals and alloys through accurate diffusivity measurements over a wide range of 
temperatures, including the proximity of the materials melting points. Specifically, it was driven 
towards developing a methodology (and subsequent flight hardware) to enable several diffusion 
coefficient measurements (i.e., at several different temperatures) to be performed using a single 
sample. 

The Liquid Metal Diffusion (LMD) was funded as a Flight Definition Project in February 1993 
in response to NRA 91-OSSA-20 (Microgravity Science and Applications Division). The Science 
Concept Review for LMD was held during April 1994. In January 1995 we were informed that 
we had failed this review and the project was change to ground-based activities only. A new 
proposal was submitted for the next NRA addressing the panels concerns. 

As part of NASA’s Risk Mitigation program, a scaled-down version of the hardware was 
funded in July of 1995 for a flight opportunity utilizing experiment on the Microgravity Isolation 
Mount. This experiment was to determine the self-diffusivity of indium at 185°C. The LMD was 
transferred to the Mir Space Station in STS-81 and returned on STS-84 (January - May 1997). 
Three, out of five, self-diffusion data sets were returned. A description of this 
experiment/hardware is included below. 

This summary is only intended to give the reader an overview of the results obtained for the 
tasks outlined in the original proposal. Research that was not published is explained in more 
detail. At the end of this report is a list of refereed publications and invited talks that were given as 
a result of this work. The reader is directed to these for further details. 

Background 

The diffusion of species in liquids governs numerous materials preparation processes of 
technological importance. In particular, a detailed understanding of diffusion in liquid metals and 
alloys is pivotal for the interpretation of various phenomena encountered in metallurgical and 
semiconductor manufacturing processes. Transport experiments in gases and crystalline solids can 
readily be designed to yield diffusion coefficients with great accuracy at normal gravity [1]. In 
liquids, however, such measurements are hampered by uncontrollable and significant contributions 
from convection. 

Simple scaling illustrates the difficulty of obtaining purely diffusive transport in liquids: In a 
system with a diffusivity of 10" 5 cm 2 /s and a characteristic diffusion distance of about 1 cm, the 
diffusion velocity is of the order of 10' 5 cm/s. If true diffusion is to be observed, convective 
contributions to mass transport must be small compared to the diffusion fluxes. Therefore, 
convection velocities parallel to the concentration gradient must be of order 10" 7 cm/s (10 A/sec !) 
or less. Hence, the results of diffusion experiments in liquids tend to be contaminated by 
buoyancy-driven flows. 

As a consequence, the standard deviations of liquid diffusivity data obtained on Earth range 
typically from several 10% to several 100% [2]. This high uncertainty is particularly deplorable 
since, due to the complex structures of liquids and, thus, the lack of realistic structure models, our 
theoretical understanding of diffusion in liquids is very limited. There are currently various 
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theoretical models for the temperature dependence of diffusivities in liquids. The predictions range 
from a pure Arrhenius-type behavior (exponential increase with temperature) to combinations of 
Arrhenius laws with power laws and pure power laws [2,3]- However, due to the large scatter of 
existing liquid diffusivity data, it is currently not possible to adequately evaluate these theories. 
Thus, accurate liquid diffusion data are needed to provide guidance in the development of both 
fundamental diffusion theory and transport models essential for materials processing. 

Efforts to exclude convective transport contributions in diffusion experiments by using narrow 
capillaries have revealed contamination of the results by (poorly understood) wall effects [4, and 
references in 2]. It appears as if, in analogy to surface and interface diffusion in solids, depending 
on the materials involved, lower or higher diffusion fluxes exist in a thin layer at the container 
wall. This effect sets lower limits to the useful capillary diameter. Thus, the drastic reduction of 
buoyancy-driven convection under low gravity conditions provides a promising means of studying 
true diffusion in liquids. 

Objectives and Tasks 

Within this framework the following technological and scientific objectives were proposed: 

• Development of a technique for dynamic in-situ measurements of diffusivities and their 
. temperature dependence in melts, with much higher efficiency than current (one sample/one 

data point) approaches permit. 

• Development of a flight-certified hardware package (GAS container) to perform such 
diffusivity measurements automatically, with relatively high flight frequency. 

• Establishment of a large definitive database for the temperature dependence of self-diffusivities 
in liquid metals and alloys, in order to further the development of the theory of diffusion in 
liquids and our understanding of numerous diffusion processes underlying materials 
processing. 

• Exploration of the possibility to approximate low gravity diffusion conditions in conducting 
liquids on Earth through the application of magnetic fields. 

WORK PERFORMED UNDER THIS CONTRACT 

Measurement Methodology (task 1) 

The experimental methods for diffusivity measurements can be grouped into direct or tracer 
techniques, and indirect or probe techniques [3,5,6]. The direct techniques employ radioactive or 
stable isotopes and diffusivities are deduced from Fick’s laws. In the indirect methods, some 
relaxation parameter associated with the macroscopic transport process is determined, from which 
the macroscopic diffusivities are deduced through model considerations. Due to this limitation we 
chose to apply a direct technique for the measurement of self-diffusivities. 

There are several direct techniques for the determination of diffusivities. They include the 
capillary-reservoir, long capillary or diffusion couple, the shear cell and the diaphragm cell 
techniques. In all of these techniques, diffusion coefficients are typically obtained by heating the 
diffusion couple through the melting point to the specified temperature, holding the sample for a 
define period of time, and then quenching the sample. The diffusant concentration is determined 
throughout the sample, and the diffusivity calculated from these concentration versus distance 
values. 
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This approach is both inefficient, yielding only one data point per sample, and prone to 
additional convective contaminations [see below] of the data due to the melting, and subsequent, 
freezing of the sample. Therefore, a primary goal of this project was to develop an in-situ, real- 
time, i.e., a radiotracer, measurement technique that would allow several diffusivity points to be 
determined in a single sample. Note that this requires being able to calculate the diffusivity from an 
arbitraiy starting profile, not only a “stepped” starting profile. 

Early methodology investigations 

We planned to employ a more “traditional” algorithm using 15-30 detectors along the length of 
the sample to monitor the concentration distribution. The initially solid, cylindrical sample contains 
a radioactive isotope at one end. After melting, radiation escaping through small bores in an 
isothermal liner/radiation shield is monitored via a chain of detectors. The diffusivities would then 
be determined from the temporal records of evolving concentration profiles through measurements 
of radioactive tracer emission. Diffusivity data could be gathered over a range of temperatures in a 
single experiment. 

Utilizing the different radiation absorption behavior of different photon energies, we could 
also investigate the significance of the “wall” effect. This effect is believed to contaminate 
diffusion studies in narrow capillaries used to suppress convection at normal gravity. 

In this approach, no radiation is received from the very ends of the diffusion sample. In 
addition, several consecutive profiles are determined at a given temperature. Then the temperature 
is changed, and the recording of consecutive profiles resumed; etc. Hence, an algorithm was 
needed that does not require data from the diffusion cell boundaries and is not limited to specialized 
initial conditions (such as a step function): 

Starting from an initial arbitrary concentration profile C(z), Fick’s second law 

C z (z,t) = DC a (z,t) (1) 

is solved for the diffusivity, in the interval a— »b such that 0 < a < b < L, with L the length of the 
diffusion cell. 

Define an arbitrary function f(z) that meets the following restrictions: 

. f I (a)=f I (b) = 0 (2) 

Thus, by integrating over n consecutive concentration profiles (using a third order spline to the 
discrete data) D is obtained for any initial conditions and without use of the physical boundary 
conditions. 

Extensive numerical modelling was performed using a variety of f(z)’s and concentration 
profiles. When using ideal (temporal) concentration profiles the input diffusivity was always 
recovered within a few 0.1%. However, on applying more realistic uncertainties (1-3 standard 
deviations) to the concentration data the uncertainty in the resulting diffusivity was greater than 
10%. These results were independent of the number of detectors used (10-30) 

Therefore, this “traditional” approach was abandoned and a two detector approach initially 
used by Codastefano, Zanza, and DiRusso [7] was tested. 
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Codastefano methodology 

Codastefano et al., showed that the higher order terms in the Fourier expansion of the 
diffusion equation canceled and the remaining one were small (compared to the First two) when the 
two detectors were maintained at L/6 from each end of the sample, i.e., 

i\ = L/6andz 2 = 5L/6. (3) 

see fig la. Then, under these conditions, the diffusivity D can be determined “directly” from the 
slope of the count difference versus time, 

ln[«i(r) - «2(0] = constant - (it/L) 2 D t , (4) 


where L is the sample length at the measurement temperature n,(t) and n 2 (t) are the total counts 
(over a given time) at measurement time t. The intensities n\ and n2 are assumed to be proportional 
to the concentrations C\ and C 2 , respectively. The characteristic shape of the signal traces «i(f) 
and 02 (f) associated with the spreading of the diffusant is plotted in Fig. lb. Obviously, before the 
first measurement, Co(z) must have spread enough to provide a significant signal at both detector 
locations. 

The concept underlying our setup was originally developed for diffusivity measurements with 
gaseous krypton. In these experiments the overall cell length was fixed; thus maintaining a single 
L. Measuring diffusivities in liquid samples over a (wide) temperature range with a single sample 
the overall sample length L, of course, changes throughout the experiment. Thus, for our 
application of this technique, additional detectors are required to maintain the advantageous L/6 and 
5176 measurement locations However, we found two pairs of collimation bores and detectors are 
sufficient. They are positioned such that, taking into account the thermal expansion of the sample, 
support structure and radiation shield, one pair fulfills condition (1) at the lowest measurement 
temperature and the other pair at the highest. As a consequence, for all intermediate temperatures, 
there is always a collimation bore above and below the exact positions required. (Note that the size 
and fixed position of the detectors are chosen such that most of the radiation emanating from a 
“moving” collimator remains detected). 

For the corrections we assume that A and B, respectively, are the diffusant concentrations at 
the collimator positions z and z + which bracket the concentration C at the accurate measurement 
position z m defined by Eq. 3. Then, noting that, at the stage where sufficient signal is obtained at 
the second collimators, the curvature of the concentration profile is small, we approximate the C at 
jt m by the linear interpolation 


C = 



(5) 
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Fig. 1. (a) Schematic presentation of the evolution of the concentration profile and location of the 
measurement locations, (b) Timetraces of the signals at the two measurements locations, and 
corresponding presentation of the signal difference according to Eq. (2). 


Extensive numerical computations showed that this methodology is very robust. The input 
diffusivity could be recovered to with a few ) 0.1% even with ± 2 standard deviations on the input 
concentration (count rate) data. 


Hardware Design and Construction (task 2) 

Choice of radiation detectors 

Since the experimental hardware developed for this project would be flown in a spacecraft the 
initial requirements for the radiation detectors were: 

• operational with cooling only to room temperature, 
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• a sufficiently high absorption coefficient, for the photons of interest, such that no more than 
1 cm of material would be necessary, 

• spectral resolution of better than 5 keV, 

• low energy pedestal of no more than 15-20 keV 

These requirements “limited” us to using either solid state detectors (i.e., CdTe, Hgl 2 , Si) or 
inorganic scintillators such as Tl-doped Nal or Csl. Silicon was rejected because of its low 
absorption coefficient and low spectral resolution. We therefore tested the following materials, 
high-pressure bridgman-CdTe, Cl-doped CdTe, CdZnTe, high-pressure bridgman CdZnTe, Hgl 2 , 
and Tl-doped Csl. In order to test for the full range of photon energies (both high and low) that 
we could use in our self-diffusion studies these materials were tested with I 129 (30, 40 and 65 keV 
photon energy) Cd 109 (22 and 88 keV), Am 241 (60 keV) and Se 75 (121, 136, 264, 280, and 400 
keV). The selection criteria between the various detectors were the peak to valley ratios, the full- 
width at half-max (FWHM) to the peak height, location of the low-energy pedestal, and the 
spectral resolution. These materials were tested at room temperature up to 65 °C. 

Based on the above selection criteria, CdZnTe was selected as the detector material. As a 
supplier we chose eV Products, Inc. (Saxonburg, PA) which had done (and continues to do) 
extensive work in miniaturizing the detector circuitry and noise reduction. 

Diffusion sample and ampoule 

The considerations for the optimum length of the diffusion couple for a Codastefano based 
methodology is different than for long capillary technique cells. In the Codastefano technique the 
diffusivity is determined from the signal difference between the two (L/6 and 5L/6) detector 
locations. Also, before valid measurements can begin, both the initial maximum “wave” of 
diffusant must pass the first detector and a (statically) significant signal level must be present at the 
second, far, detector. Thus, there is an optimal sample, and diffusant, length, such that 
measurements are valid shortly after the maximum wave has passed the first detector and such that 
the diffusant concentration is equalibrated (fully diffused) over a long enough time to allow for 
multiple determinations. In addition, the sample diameter is optimized to allow for a reasonable 
(radioisotope) activation parameter, i.e., if the isotopically enriched portion is too small, a much 
long activation time is required. Finally, the sample diameter should be large enough to allow us to 
investigate the “wall effect”. Following extensive numerical modelling we found that a 3 mm <]) 
sample with a 29 mm long native abundance section and a 1 mm radioactive tracer (diffusant) 
portion sufficiently fulfilled the criteria. 

Then, as indicated in Fig. la, we use an initially solid cylindrical diffusion sample which 
consists mostly of inactive material (solvent) and a short section of activated isotope, as the 
diffusant, located at one end. After melting of the sample and heating to a uniform measurement 
temperature T, the evolution of the diffusant concentration distribution C(z) is monitored through 
the intensity of the radiation received through two bores (collimators) in a radiation shield. 

BN Ampoule design 

The BN ampoule consists of four pieces and a graphite spring, see figure 2. The four pieces 
are the ampoule body, bottom plug, plunger and screwcap, which are constructed from 
Carborundum grade AX05. This high-purity grade BN is not wet by any of the metals under 
consideration for this work. The ampoule body I.D. is reamed to 3.00, with tolerances of -0 and + 
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0.01 mm to minimize errors in the (total) sample length. The ampoule body O.D. varies depending 
on the wall thickness of the particular outer cartridge. The bottom plug is precision machined to 
3.05 ± 0.02 (J) and 5.0 ± 0.02 mm long. It is press fitted into the body, providing "square" corners 
at the sample bottom. The plunger is individually fit to less than 0.02 mm clearance to the ampoule 
body. Close tolerance between the plunger and the ampoule bore prevents liquid metal from 
escaping. The low-force graphite spring pressing against the plunger 



Fig. 2. Cross-section of diffusion sample cartridge. 

ensures a flat interface at both ends of the sample. This ampoule is vacuum-baked at high 
temperature to remove gaseous contaminants. 

An outer cartridge provides a vacuum or gas seal to minimize oxidation or transport of the 
diffusion couple. The specific material used varies depending on the temperature range and 
radiotracer employed. Up to the present time, silica glass, aluminum and stainless steel have been 
used for these cartridges. 

Sample preparation 

' The native abundance section is prepared from 6N starting material. It is vacuum-cast (10~6 
mm Hg) in a 2.9 mm diameter BN mold to remove residual oxides, voids and gaseous 
contaminates. It is then removed, trimmed and weighted to 0.1 mg. This section is then loaded 
into the diffusion ampoule itself and vacuum cast again to approximately 100 °C above the planned 
processing temperature. The final vacuum casting minimizes any voids and residual oxide in the 
material. 
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The isotopically enriched (diffusant) section is prepared into a flat disk using a boron nitride 
ampoule consisting of three parts machined to tight tolerances. After loading and assembling the 
ampoule the isotopically enriched material is vacuum melted, again to approximately 100 °C above 
the planned experiment temperature. It is then removed under inert gas, weighted to 0.1 mg and 
sealed in a silica glass tube for neutron activation. The activated section is removed from the silica 
tube, under inert gas, and loaded into the ampoule. After replacing the plunger and spring the BN 
ampoule is mechanically closed with the screw lid and inserted into the (square-bottomed) outer 
cartridge. 

Isothermal liner / radiation shield 

The Isothermal liner/radiation shield was to function as both a zone leveler, to minimize any 
temperature nonuniformities across and along the diffusion sample, and as a collimator for the 
emitted photons. As originally conceived LMD was to be flown as an “shuttle-independent GAS 
cannister experiment. Thus, both experiment mass and the power required for operation were to 
be minimized. Gold was chosen as the liner/shield material due to its high thermal conductivity, 
high radiation absorption and ready machinability. The initial liner/shield design yielded a cylinder 
approximately 21 mm diameter (7 mm wall) and 90 mm long. Sixteen 0.67 mm collimators were 
bored in a spiral pattern along the liner/shield axis. In order to minimize heat transfer, and 
therefore, the input power required, considerable effort was devoted to reducing the convective and 
radiative heat transfer. As presented at the LMD Science Concept Review (SCR) the liner/shield 
was suspended by a wire “spokes”. A series of gold-flashed thermal shields minimized the 
radiative losses. With this design we were able to maintain sample temperatures of 800 °C with 20 
watts. Assuming 28 V power (or batteries) this represented a current draw of less than 1 A. 

However, in the first diffusion measurements with radiotracers excessive background through 
the shield overwhelmed the signal (radiation) from the collimators. Thus, it was obvious that 
either the shield wall thickness was not sufficient or that the collimation hole diameter was too 
small; or both. 

A series of experiments was run to determine the optimum collimator geometry. Both the 
collimation bore diameter and the shield thickness were varied, between 1 and 3 mm, and 13 to 39 
mm of Pb, respectively. A small collimation bore diameter gives high spatial resolution but a low 
signal-to-background level; for a large collimation hole the converse is true. A thin shield 
thickness (low bore length) has minimum scatter in the bore but high direct radiation penetration 
and fluorescence. The results of these tests indicated that a 2 mm diameter collimation hole and a 
16 mm bore length would provide the optimum signal-to-background level for diffusion 
measurements. Thus, the gold liner/shield that has been used for all measurements to date under 
this contract has been 41 mm diameter, with an approximately 9 mm center bore (see below), 70 
mm long. 2 mm diameter collimation holes are bored through the liner/shield at the required L/6 
and 5L/6 positions. 

Choice of Radioisotopes to be Used (task 3) 

As stated above, an original objective of the LMD project was to development a measurement 
methodology and flight hardware. Beyond suggesting individual elements or alloys for 
investigation little justification was given for specific choices. The Science Concept Review panel 
found this to be scientifically unsound. In response to this criticism, in preparation for submitting 
to the next available NRA, we systematically searched the literature for possible liquid structure 
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features that may be responsible for specific diffusion mechanisms. From this we selected a series 
of elements for study based on their “class-like” self-diffusion behavior. 

Class-like self-diffusion behavior of molten elements 

In spite of the large uncertainties in diffusivity data, enough insight has been obtained, in 
particular through a few measurements spanning wide temperature ranges, to recognize some 
class-like behavior that can be associated with temperature -dependent structure features of liquids. 
D(T) data are conveniendy tabulated in terms of the two parameters Do and Q in the relation 
D = D 0 exp(-Q/RT), where Q is understood only as an apparent “activation energy”, without the 
intent to support activated state theories. Such compilations of data, and corresponding plots of 
In D vs. 1/T readily reveal a few cases of clear “non- Arrhenius" behavior, i.e. discontinuous, steep 
increases of the diffusivity with temperature. 

Tin, for instance, shows a sharp increase of D with rising T around 1030 °C [8]. This has 
been associated with the dissociation of a second structure in the liquid [8]; see below. This 
dissociation is indicated by the disappearance of a subsidiary peak in plots of the scattered X-ray 
intensity versus scattering vector [9] and the appearance of a peak in the differential thermal 
analysis curve [10]. Similarly, in liquid gallium, below 70 °C the Arrhenius plot yields a diffusion 
“activation energy” of 1.1 kcal/mol, while at higher temperatures it is about twice this value [11]. 
Correspondingly, a subsidiary X-ray peak disappears at around 100 °C [12]. Mercury shows a 
similar behavior, with the disappearance of the second liquid structure and the sudden increase in 
D(T) occurring at around 0°C [10]. Liquid tellurium, the only elemental liquid semiconductor (see 
also [13]), on the other hand, possesses a sharp decrease in the slope of D(T) with increasing 
temperature around 500 °C, i.e. 50 °C above the melting point [14], Again, neutron scattering data 
suggest a structural transition around this temperature [15]. The initial steep slope in D(T) has 
been interpreted in terms of the dissociation of chain-like structures and, thus, increased thermal 
mobility of molecules in liquid tellurium, that is essentially completed around 500 °C [14], In 
contrast to these complex systems, the alkali metals as well as copper and silver show rather 
monotonic D(T) behavior [3]. 

The above cases suggest that for liquid elements for which detailed diffraction and, thus, 
structure data exist, valuable predictions of the D(T) can be made. Such structure data exist for a 
large number of liquid metals and semimetals [16-18]. Most simple liquids and liquid metals 
exhibit a sharp main peak in plots of scattered intensity versus scattering vector. Liquid metals 
exhibiting a symmetrical sharp maximum (Saweda type-a) include the alkali metals. Mg through 
Ba, copper, silver, gold and the transition metals. Symmetrical main peaks, which are also 
obtained from hard-sphere-interaction diffraction models [18], indicate complete randomness in 
atom arrangement, i.e. no preferred configuration of atoms, however transient [16]. 

In several liquid metals, the first scattering peak is asymmetrical (Saweda type-b [18]), and in 
some cases has a pronounced shoulder on the high-angle (lower atomic distance) side of the peak 
(Saweda type-c). This group includes, in addition to Sn, Ga and Hg mentioned above, Zn, Cd, 
Si, Ge, Sb, and Bi, with possible minor asymmetries in the peaks of In [19] and Tl [20]. Note 
that both Ge and Si, which are semiconductors in the solid state, turn metallic on melting [13]. In 
the non-metallic liquids Se and Te, which show predominantly covalent (homopolar) bonding, this 
subsidiary peak is particularly pronounced. Thus the probability of finding an asymmetric 
maximum is higher in elements of the higher groups, in which the solid-state bonding is 
predominantly non-metallic [21], The distance values corresponding to these subsidiary peaks 
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often represents the nearest-neighbor distance characteristic of the covalent bond in the solid state, 
e.g. gray tin in the Sn-system, while the main peak closely correspond to the distance in the 
metallically-bonded solid, e.g. white tin. This suggests that in these liquids there is a short-term 
localization of valence electrons in bound states between pairs or groups of neighboring atoms 
within an otherwise completely random metallically bonded matrix. With increasing temperature, 
the formation of short-lived clusters becomes less likely, leading to the disappearance of the 
subsidiary diffraction peaks. Even more pronounced than in Se and Te, liquid sulfur (an insulator) 
shows on heating several structural transitions, between various molecular ring and chain 
configurations, over a range of more than 100 °C above the melting point [22,10] 

High/Low energy photons for “ wall effect” investigations 

In addition to its real-time feature, this isotopic labeling technique permits investigations of the 
perceived “wall effect” mentioned above. By employing an isotope which emits photons at two 
sufficiently different energies and, thus, different self-absorption behavior, transport in the bulk of 
the sample and near the container wall can be distinguished to some extent. For example, based on 
attenuation data, for the 24 keV and 190 keV photons of 114m In [23] we have calculated the 
fraction of the total radiation received at the two different emission energies outside a 3 mm thick 
sample vs. the distance of the emitter from the surface of the sample; for details of the one- 
dimensional slab model used see Refs.[24], Note that the 24 keV photons received are predicted to 
originate in essence only from a 300 |im deep surface layer. The 190 keV photons, on the other 
hand, stem from throughout the whole sample slab. This behavior was experimentally confirmed, 
for details see Jalbert, Banish, and Rosenberger at the end of this report. Thus, in selecting 
elements for self-diffusion studies consideration was given to isotopes that emitted (multiple) 
photons both high enough in energy to be representative of “bulk” diffusion and low enough to 
represent transport in the “wall” region of the sample. 

Systems to be used in diffusivity measurements 

Based on the above considerations we proposed investigating the following systems. 

• Simple liquids: 

Indium which may show some weak nonmonotonic behavior, has a combination of low and 
high energy photons that can be utilized for quantifying wall effects, a conveniently low 
melting point and a wide liquid temperature range. 

Rubidium has been studied over approximately 80% of its liquid state. However, the 
temperature dependence of Rb has been reported as both T^ and Arrhenius. 

Aluminum based on neutron diffraction data aluminum is expected to have a monotonic 
temperature dependence. However, aluminum self-diffusivity values are not available. 

• Complex liquids: 

Cadmium has two advantageous photons for studying bulk and wall transport. 

Gallium with low melting point and wide liquid range with x-ray diffraction data showing a 
pronounced asymmetry in the first peak. 

Tin which has been extensively investigated on Earth and to some extent in space [2,8,9], 
producing controversial results. It also has two particularly advantageous photons for 
studying bulk and wall transport. 
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Tellurium and selenium, which should provide interesting features in D(T), due to the 

multiple structures that are known to exist in their liquid states. Te has the advantage of 
two distinctly different energy photons and Se that of a lower melting point. Note that 
the energies of Te also appear well suited for wall effect investigations. 

Convective Effects on Diffusivity Measurements (task 4) 

Numerical modelling of magnetic fields to suppress convection 

The scaling considerations presented above show that convection velocities in typical self- 
diffusion experiments must be reduced to much less than 10’ 5 cm" 1 . In principle, this can be 
achieved in liquid metals using a steady magnetic field to damp the flow. In order to assess the 
magnitude of the field required, we have computed buoyancy-driven flows in a model diffusion 
cell geometry with and without applied magnetic field. 

The 2D Navier-Stokes-Boussinesq equations were solved for flow driven by small lateral 
temperature gradients in a vertical cavity of 5 cm length with a width of 1 and 3 mm, respectively. 
As illustrated in Fig. 1, two types of lateral temperature boundary conditions were considered. In 
case 1, a uniform temperature difference of either 0.1 or 0.01 K was applied across the slot. In 
case 2, a “cold spot” of about 0. 1 K was assumed at the upper right comer. Transverse and vertical 
magnetic fields were considered and it was assumed that generation of secondary magnetic fields 
by the flow could be neglected. The two-dimensionality of the problem allowed us to neglect the 
electric field [25]. 

The steady Navier-Stokes Boussinesq equations governing the transport of heat, mass and 
momentum are written in dimensionless form as 

U ^ii = - Vp + PrV 2 u - RaPr0k - Ha2 (uxB)xB, (6) 

ot 

Vu =0 ' (7) 

and 

u- ve = v 2 e, * (8) 

dt 

* ■ t 

where u(x,t) , B and 0(x,t) are the dimensionless velocity, magnetic field and temperature, 
respectively, and Ra = pgATd^/KV, Pr = v/k and Ha = Bdfa/p.) 1 / 2 . The physical properties P, 
AT, d, k, c, (i and v are, respectively, the thermal expansion coefficient, the maximum deviation 
from the desired uniform temperature, the width of the slot, the thermal diffusivity, the electric 
conductivity, and the dynamic and the kinematic viscosity. Ra, Pr and Ha are the Rayleigh, 
Prandtl and Hartmann numbers, respectively. 

The governing equations were recast in stream function vorticity form and discretized using a 
Chebyshev pseudospectral collocation method with an Adams-Bashforth/second backward Euler 
temporal scheme [26,27]. The discretized problem was solved using an influence matrix method 
[26-28]. It is well known that at high values of the Ha the flow is confined to narrow boundary 
layers with thickness proportional to Ha -1 / 2 for walls perpendicular to the field lines and Ha -1 for 
walls parallel to the field lines [25]. Thus, we needed to ensure sufficient spatial resolution to 
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Case 1 Case 2 



Figure. 3. Sketch of problem domain for both case 1 (lateral temperature gradient extends 
throughout the domain) and case 2 (lateral temperature gradient on the right-hand wall only 
extends over 2 cm). 

capture the boundary layers. Collocated Chebyshev polynomials' were distributed such that for 
N+l points in a given coordinate direction xj = cos(iic/N), i = 0,1..., N.. We used 33x85 
collocation points. This gave us 6 points inside 1.225 mm for the horizontal boundary layers (for 
the 3 mm cell width) and 6 points inside 0.477 mm for the vertical boundary layers. Given that the 
points are more closely spaced as one approaches the boundary, this appears to provide sufficient 
resolution for these flows. 

Using the physical properties of liquid tin, we obtained solutions for both temperature 
boundary conditions. The lateral magnetic fields were found to be the most effective for damping 
the flow. This is because the main damping effect is on the velocity component perpendicular to 
the field lines. In the absence of a magnetic field, the maximum velocity is the vertical component, 
w. A vertical field is not as effective in reducing the w-velocity since it only directly damps the u- 
component. Nevertheless, in the bulk of the domain, flow velocities are diminished and the flow 
is squeezed close to the boundaries. Three sets of results for lateral fields are presented in Table 1 . 
For case 1, both cell widths were considered. Since for case 2, the 3 mm width velocities at high 
Hartmann numbers were similar in magnitude to those in case 1 we did not carry out calculations 
for 1 mm widths for this case. The flow velocities obtained in the absence of a magnetic field are 

1 The Chebyshev polynomials are defined as T n (x) = cos(n arccos x), n = 0,1,..., and -1 < x < 1. 
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about the same as those predicted by the approximation of Batchelor [29], i.e. ~ 10' 2 RaK/d. For 
Ha = 0 the calculated results scale linearly with a decrease in temperature or with gravity level. 

Summary 

We have calculated the convective flow in diffusion cells arising from small horizontal 
temperature gradients with and without steady, transverse magnetic fields. Only for high fields, 
on the order of 1 Tesla, could velocities be brought down to the order of 10" 5 cm s' 1 or less. At 
high lateral fields the maximum velocity is the horizontal component, which is typically twice the 
vertical maximum. Since typical diffusive speeds will be on the order of 10" 5 cm s*l (based on v~ 
D/L, where D is the diffusivity and L= 5 cm), it appears that, even for lateral temperature 
nonuniformities as low as 0.01 K, unless a combination of thin capillaries (less than 1 mm width) 
and large magnetic fields (in excess of 1 Tesla) are used, the maximum flow speeds cannot be 
reduced significantly below diffusive speeds. 


Table 1: Convection speeds 

in diffusion cell model with and without transverse magnetic field. 

Case 1 


■ Sill 

Maximum velocity (3 nun) : 

!> B 

Maxirnum velocity (1 Win»); : 



m 

[cm s* 1 ] 

■Bill 

• [cm s' 1 ] J= :: 

AT = 0.1K 

Ha 


u 

w 


u 

w 


0 

0 

1.13xl0- 2 

1.87x10-2 

0 
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2.08x10-3 
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3.26xl0- 3 

6.57x10-3 

0.29 
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50 

0.5 

6.5xl0’ 4 

1.3x10-4 

1.47 

3.32x10-5 
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100 

i 

1.29xl0- 4 

5.01x10-5 

2.95 

1.25x10-5 

6.09x10-6 
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1.47 

7.60xl0- 5 
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0 
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10 
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3.52x10-5 
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50 
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1.24x10-6 
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Numerical modelling of buoyancy -driven convection in diffusion cells 

To provide more quantitative guidance for experimental work, we numerically modeled 
diffusive-convective transport in vertical cylinders. The dimensions of the cylinders were chosen 
to correspond to typical diffusion experiment geometries. Thermal boundary conditions with 
slight asymmetries were imposed to simulate small deviations from isothermality. Since the 
objective is to simulate convective effects on self-diffusion, the diffusant and sample were assumed 
to have the same mass density. The three-dimensional time-dependent transport of the diffusant 
(initially located at one end of the sample) was simulated for the physical properties of liquid 
indium. 

The liquid phase was assumed to be a Boussinesq fluid contained in a closed vertical circular 
cylinder of length L and diameter 2R. Gravity acts along the cylinder axis. The temperature at the 
vertical wall of the cylinder has the form: T(r,0,z) = (T, -T 2 )/2 -(T, + T^cosS)/^ + AT z z/L. Here 
T, - T 2 is the maximum lateral temperature difference across the cylinder at any height z, and AT Z is 
the temperature difference between the bottom and the top. The temperature of the latter are T(r, 0 , 
0) = (T, -T 2 )/2 - (T1 + T2) r (cos0)/2R and T(r, 0, L) = T(r, 0, 0) + AT Z , respectively. The 
governing equations for the time-dependent transport of heat and momentum in the system were 
solved using finite difference on an equally spaced mesh and on a variably spaced mesh using the 
code CFD2000 which is based on the finite volume method. 

Summary 

The simulation yielded various results of practical consequence. We found that, even for 
temperature nonuniformities Tj -T 2 as low as 0.1 °C, convective transport in capillaries of 1 mm 
diameter can largely exceed diffusion. Furthermore, the convective contribution to transport may 
be increased by “stabilizing” AT z ’s. Hence, the widely used approach of keeping the top of 
diffusion samples slightly warmer than the bottom appears to be a dubious means of “preventing 
convection”. Most striking, we found that diffusive-convective transport can result in profiles of 
the radially averaged concentration values that are analytically indistinguishable from pure diffusion 
profiles. This finding casts considerable doubt on various geometrical criteria previously used to 
demonstrate the absence of significant convective contamination diffusivity data plots. In 
summary, reduced gravity conditions are essential for definitive diffusivity measurements in liquid 
phases of macroscopic dimensions. For details of this work see (Alexander, Ramus and 
Rosenberger). 
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Risk Mitigation Flight Opportunity 

A simpler version Liquid Metal Diffusion experiment was selected for a flight opportunity 
under the NASA/Mir Glovebox/MIM risk mitigation program. The scientific and technological 
goals of this phase were: 

• Determination of the effects of steady residual gravity and g-jitter aboard spacecraft on 
liquid diffusivity measurements. 

• Determination of the influence of g-level on possible “wall effects”. 

• Test of anchoring technology to prevent void/bubble formation in liquid metal diffusion 
samples at low gravity. 

• Tests of a sample exchange technique for multiple runs of a diffusivity measurement setup 
during a spaceflight. 

• Test of radiation detector performance with respect to space background radiation level and 
sample activity requirements. 

The LMD hardware utilized the (Canadian Space Agency) Microgravity Isolation Mount 
(MIM) as part of the Mir Space Station. The MIM hardware provided a unique opportunity for 
testing the effects of residual gravity and g-jitter on liquid diffusion sample. It provided for 
vibration isolation and defined g-inputs up to 10' 3 g 0 ; as well its’ latched configuration; where g- 
gitter and residual accelerations would be transmitted to the samples. 

To take full advantage of the MIM capabilities, five ll4m In/In metal self-diffusion samples were 
flown. All samples were to be processed at 185 °C. (Note, initially the samples were to be 
processed at 200 °C. However, because of concerns from the Russian Space Agency, the 
processing temperature was decreased to 185 °C) The operating characteristics of the MIM during 
sample processing were to be: 

• One sample with the MIM latched down, 

• two samples with isolation only, 

• one sample with 10’ 3 to 10 2 go pulses, at approximately 0.01 Hz, normal to the sample axis, 

• one sample with 10' 3 to 10' 2 g 0 pulses, 0.01 Hz, perpendicular to the sample axis. 

The design and construction of the flight hardware employed the same systems and 
technology used for ground-based experiments. Containment and electronic redundancy 
requirements that were placed on the hardware did not change its functionality. Additional samples 
were contained in a hand operated storage carousel not typically used on the ground. Figure 4 is a 
scaled cross-section of the flight hardware indicating the major components. 
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Figure. 4. Liquid metal diffusion apparatus cross-section: (1) diffusion sample in cartridge, (2) 
heated isothermal liner/radiation shield, (3) CdZnTe detector pairs, (4) preamplifiers/energy 
discriminators, (5) circuit board, (6) carousel with four additional sample cartridges, (7) lead pig, 
(8) sealed aluminum housing. 

Application of technique to indium: ground-based and low-gravity results 
Indium was chosen as the first element for investigation due to its: 

• low melting point of 1 56 °C, 

• high boiling point of approximately 2000 °C, thus, giving a very wide liquid range, 

• radioisotope having low (24 keV) and high (190 keV) photons with a high enough 
branching ratio such that activation was relatively straight forward. The energy range of these 
photons was such that wall effects could be investigated. 

Diffusivity determinations were performed at a single temperature as a verification of the 
method and hardware, both on Earth and aboard the space station Mir. In the ground-based tests 
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we also investigated the effect of the initial source activity levels on the signal-to-noise ratio. All 
experiments were run until diffusant/emitter uniformity was obtained throughout the sample. 
Typically this requires 100 hours at the low temperatures used in these experiments. 

The first five ground-based runs were performed at temperatures near 200 °C. In these early 
experiments the temperature control with a manually adjusted DC power supply was not optimal. 
Indium samples with room temperature lengths of 30 mm and different initial source thicknesses 
and activities were used. Radiation counts were taken continuously and summed every 10 
seconds. Diffusant source parameters and apparent diffusivity results are summarized in Table 2. 


Table 2. Summary of ground-based results at 200 °C. 



Self-diffusivities in indium [10"^ cr 

Left detectors 

190 keV 24 keV 

i 2 /s] 

Right detectors 
190 keV 

24 keV 

53 

2.44±0.16 

(3.74)* 

(only one pair used) 

75 

2.60±0.17 

2.4710.16 

2.6410.17 

2.4310.16 

150 

2.46±0.16 

(3.69)* 

2.4610.16 

(2.36)* 

600 

2.14±0.14 

2.1610.14 

2.2110.14 

2.2410.14 

5000 

2.3310.15 

2.0410.13 

2.2410.14 

2.0610.13 


* ground loop/stray capacitance problems: calibration shifts during experiment. 


After a thorough examination of the raw data and calculated diffusivities an initial activity of 1 
mCi was chosen as sufficient for single temperature measurements. Using this activity, samples 
were prepared for the liquid element diffusion (LMD) experiments on space station Mu- 

Four of the five diffusion samples delivered to Mir were run. Data from three tests were 
returned. Data from the second sample was lost due to software problems; the fifth sample was 
not run due to Mir time constraints, i.e., defined input g-pulse runs were not carried out. 
Therefore, instead of the MIM operating conditions stated above two samples were run in latched 
mode, i.e., without g-jitter isolation, samples 1 and 3. Run 2 was performed with active isolation. 
Table 1 shows our low temperature (approx. 200 °C) indium self-diffusivity values both at normal 
and low-gravity (less than lO' 4 - 10' 6 g 0 ). 
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Table 3. LMD Mir results. 


Run 

Temperature 

[°C] 

Self-diffusivities in indium [] 

Left detectors 

190 keV 24 keV 

0-5 cm^/s] 

Right detectors 
190 keV 24 keV 

1 (latched) 

2(isolated) 

3(latched) 

186.5 ± 0.8 
185.0 + 0.4 
186.5 1 0.4 

1.9910.13 2.0710.13 

2.0910.13 1.9810.13 

2.0810.13 2.0910.13 

1.9810.13 2.0610.13 

2.0710.13 * 

2.0510.13 * 


* ground loop/stray capacitance problems: calibration shifts during experiment. 


Discussion of low temperature results 

The average (normal gravity) values of the self-diffusivity of In 114m /In at 185 °C are; 

• 2.22 ± 0.159 (for 190 keV photons) and, 

2.11 1 0-154 x 10' 5 cm 2 /sec, (for 24 keV photons), 

(The normal gravity values in Table 21 were corrected to 185 °C by the factor 1 x 10 7 cm 2 /sec/K.) 
The average, low-gravity, indium self-diffusivity values are; 

• 2.04 1 0.047 (for 190 keV photons) and, 

• 2.04 ± 0.044 x 10' 5 cm 2 /sec (24 keVphotons). 

These low temperature experiments established that this technique and hardware will provide 
accurate diffusivity data in element liquids. Several aspects are noteworthy. First, within the 
experimental uncertainty, the results are reproducible. Furthermore, the agreement between the left 
and right detector pairs, see Tables 2 and 3, indicates that at these low temperatures either 
convective transport contributions are insignificant or, less likely, affect the concentration 
distribution in the measurement plane rather symmetrically. The latter could result from a 
convection roll that is oriented normal to the plane of measurement. The reduced spread of the 
low-gravity data, while the D,’s themselves are not that different, are difficult to unambiguously 
explain. Similar results have been noted by other authors [30]. While, there may be some 
convective contributions (in the ground based data) to the apparent diffusivity, they are probably 
nonsteady and, therefore, difficult to characterize fully, or eliminate. 

Given the low scatter and reproducibility of the reduced-gravity data we report these values 
with good confidence. Therefore, we do not see a need to measure these values again. The 
reproducibility between the isolating and latched runs indicates, coupled with our numerical 
modelling, that the LMD cell was isothermal to better than 0.1°C, which was a design goal. 

Overall, these results are on the low end of the range of published values for D/ n (200°C) 
which are from (2.2 - 3.8 ) x 10"5 cm^/s. Finally, it should be noted that for all data, transport in 
the bulk of the sample and near the container wall were indistinguishable, i.e., no ’wall effect' was 
detected. 
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Risk Mitigation Results 

The overall results from this flight opportunity were very positive and for most aspects the 
experiment was a success. In addition to the low temperature indium self-diffusivity results: 

1) the sample exchange and positioning mechanism worked properly, 

2) the radiation detectors worked well within specifications, 

3) the chosen sample activity was acceptable for the (Mir) background radiation environment, 
and 

4) the sample containment/restraint system work properly. 

Summary 

During the course of this work, most of the proposed objectives were successfully completed. 
We did not pass our Science Concept Review and, therefore, we did not build and fly hardware 
capable of measuring diffusivities of several elements over a wide temperature range. However, 
we did develop a methodology for this that has been demonstrated on the ground over a wide 
temperature range. This same Codastefano methodology was successfully used for our low 
temperature spaceflight results. Hardware components, radiation detectors, isothermal liner/shield, 
etc., were also developed during this project. Numerical modelling showing the extreme 
sensitivity of diffusion measurements to convective contaminations was completed. Finally, 
numerical modelling of the self-defeating nature of magnetic fields to suppress convective 
contaminations, in diffusivity measurements, was also completed. 
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Based on the methodology of Codastefano, Di Russo, and Zanza [Rev. Sci. Instrum. 48, 1650 (1977)], we 
have developed a technique for the in situ measurement of diffusivities in liquids at several temperatures with 
one sample. In this approach, which circumvents solidification of the diffusion sample prior to concentration 
profiling, the evolution of the concentration distribution of a radiotracer is followed in real time using two pairs 
of radiation detectors. A detailed description of the experiment setup and its optimization by computer simu- 
lations is given. In experiments with U4m In/In, apparent self-diffusivities were obtained between 300 °C and 
900 °C with an uncertainty of ±5%. By utilizing the different self-absorption characteristics of the 24- and 
190-keV photons of n4n, In, transport in the bulk of the sample and near the container wall was investigated 
independently. No difference was found. [SI 063-65 IX (98)09002-3] 


PACS number(s): 66.I0.Cb, 07.85. -m, 07.05. Fb 

I. INTRODUCTION 


Most techniques to measure diffusion in high-temperature 
liquids involve solidification of the sample before concentra- 
tion profiling [1-3]. Inevitable nonuniformities during freez- 
ing can introduce considerable concentration redistribution 
due to segregation [3]. Furthermore, in samples of low ther- 
mal conductivity, the heat of solidification can lead to sig- 
nificant temperature nonuniformities. As a result, convective 
mixing and a higher apparent diffusivity can result even in 
systems that are perfectly isothermal during the liquid diffu- 
sion process [3], 

Some radial concentration nonuniformities observed in 
samples that were solidified in diffusion capillaries have 
been assigned to “wall effects” [4]. These effects are per- 
ceived, in analogy to surface and grain boundary diffusion, 
to originate from different diffusion rates near solid walls 
and in the bulk liquid. However, as pointed out by Nachtrieb 
[3], the range of true wall effects should be limited to a few 
atomic dimensions from the wall. This is corroborated by 
analytical investigations in terms of particle density wave 
scattering by a wall [5] as well as by molecular-dynamics 
simulations [6,7]. Both approaches suggest that true wall ef- 
fects should be limited to some 50 A from the wall, i.e., to 
such a small fraction of the sample cross section that they 
should be undetectable by macroscopic diffusivity measure- 
ments. 

We have addressed the above problems and developed 
an experimental technique for real-time measurements of 
diffusivities in liquids. This radiotracer technique cir- 
cumvents solidification of the diffusion sample prior to 
concentration profiling. Furthermore, by employing an 
isotope that emits photons at (two) sufficiently different 
energies, transport in the bulk of the sample and near the 
container wall can be distinguished. In the following, we 
will present the conceptual considerations underlying this 
technique, simulations carried out towards its optimization, 
details of the experimental setup and sample preparation, 
and applications to self-diffusivity measurements in liquid 
indium. 


II. EXPERIMENT CONCEPT 
A. Methodology 

The concept underlying our setup was originally devel- 
oped by Codastefano, Di Russo, and Zanza for diffusivity 
measurements with gaseous krypton [8]. As schematically 
indicated in Fig. 1(a), we use an initially solid cylindrical 
diffusion sample that consists mostly of inactive material 
(solvent) and a short section of activated isotope, as the dif- 
fusant, located at one end. After melting of the sample and 
heating to a uniform measurement temperature 7\ the evolu- 
tion of the diffusant concentration distribution C(z) is moni- 
tored through the intensity of the radiation received through 




FIG. 1. (a) Schematic presentation of the evolution of the con- 
centration profile and location of the measurement locations, (b) 
Time traces of the signals at the two measurements locations, and 
corresponding presentation of the signal difference according to Eq. 

( 2 ). 
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FIG. 2. Emission of 24- and 190-keV photons from 3-mm-thick 
indium sample containing "“"In versus distance of the emitting 
atoms from the sample surface. 


two bores (collimators) in a radiation shield. These intensi- 
ties n | and n 2 are assumed to be proportional to the 
concentrations C, and C 2 , respectively. The characteristic 
shape of the signal traces n,(r) and n 2 (t) associated with the 
spreading of the diffusant is plotted in Fig. 1(b). 

To satisfy the requirements of the algorithm used to 
evaluate the diffusivity [8],' the radiation collimation bores 
must be positioned at 

z , = L/6, z 2 = 5L/6, (1) 

where L is the sample length at the measurement tempera- 
ture. The diffusivity D is then calculated from the difference 
of the signal traces using the relation 

ln[nj(f) — /i2(f)] = const - (ir/L) Dt, (2) 

where the constant depends on the concentration profile 
C 0 (z) at the beginning of the measurement. Since the C 0 (z) 
does not explicitly enter the D evaluation, diffusivities can 
be consecutively determined at several temperatures during 
the spreading of the concentration profile in the same 
sample. Obviously, before the first measurement, C 0 (z) must 
have spread enough to provide a significant signal at both 
detector locations. 

In addition to its real-time feature, an isotopic labeling 
technique permits investigations of the perceived wall ef- 
fect” mentioned in the Introduction. By employing an iso- 
tope that emits photons at two sufficiently different energies 
and thus different self-absorption behavior, transport in file 
bulk of the sample and near the container wall can be distin- 
guished to some extent. Figure 2 illustrates this approach. 
Based on the attenuation data for 24- and 190-keV photons 
of ll4m In [9], we have calculated the fraction of the total 
radiation received at the two different emission energies out- 
side a 3-mm-thick sample vs the distance of the emitter from 
the surface of the sample; for details of the one-dimensional 
slab model used see Refs. [10,11], Note that the 24-keV 
photons received are predicted to originate in essence only 
from a 300-/xm-deep surface layer. The 190-keV photons, on 
the other hand, stem from throughout the whole sample slab. 
Thus, using an appropriate detector circuit with energy dis- 


crimination capability (see Sec. IV C), transport near the 
wall of the sample container and in the sample bulk can be 

distinguished to some extent. 

We have experimentally confirmed the above selective 
absorption behavior. Indium disks of 100, 200, and 300 M m 
thickness were prepared from indium that had been uni- 
formly diffusion doped with 1 ,4m In. The normalized ratios of 
the 24- and 190-keV counts received from these disks agreed 
well with the predicted values; see the data points on the 
curves in Fig. 2. 


B. Collimator positions and thermal expansion 


Relation (2) is based on the specific location requirement 
for the detectors given by Eq. (1) with respect to the sample 
length L. For applications of this measurement technique at 
several temperatures with one sample, corrections are re- 
quired to account for thermally induced changes in sample 
dimension and collimator locations. For these corrections it 
is advantageous to use two pairs of collimation bores and 
detectors. They are positioned such that, taking into account 
the thermal expansion of the sample, support structure, and 
radiation shield, one pair fulfills condition (1) at 
measurement temperature and the other pair at the highes . 
As a consequence, for all intermediate temperatures, there is 
always a collimation bore above and below the exact posi- 
tions required. (Note that the size and fixed position of the 
detectors is chosen such that most of the radiation emanating 
from a “moving” collimator remains detected; see also bee. 


For the corrections we assume that A and B, respectively, 
are the diffusant concentrations at the collimator positions z 
and z + that bracket the concentration C at the accurate mea- 
surement position z m defined by Eq. (1). Then, noting that, at 
the stage where sufficient signal is obtained at the second 
collimators, the curvature of the concentration profile is 
small, we approximate the C at z m by the linear interpolation 



(3) 


ffl. EXPERIMENT SIMULATION AND OPTIMIZATION 

We have carried out detailed computer simulations to op- 
timize the experiment setup with respect to sample size and 
activity, measurement accuracy, and number of D(I ) data 
obtainable with one sample. In the following we derive the 
basic relations used; for more details on the simulations the 
reader is referred to Ref. [10]. 

A. Evolution of concentration distributions 
including sample expansion 

In a realistic simulation of the experiments, changes in 
sample length due to thermal expansion during heating to the 
first measurement temperature and heating or cooling to con- 
secutive temperatures must be taken into account. Let the 
initial concentration profile in the sample of total length L, at 
temperature T , be 

g{z) = C 0 6(z-li), 0<z<L \ . W 

with C 0 the initial concentration over the initial length of the 
activated section /, and *z-/,) a step function equal to 1 
for z </, and 0 for z>/ f . With this initial condition and the 
closed-system boundary condition 
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C.(0./) = C : (L„f) = 0, 0<f. (5) 

assuming that the diffusivity D is independent of concentra- 
tion and coordinate z, one obtains from Fick’s second law 

C,{z,t) = DC u (z.t) (6) 


a general solution for the evolving concentration distribution 
in the form [12] 


n TTZ 


C(z,/) = fl 0 +2 C0S (-^) eX P] 


n = 1 


nir 
— ) Dt 


, (7) 


with the Fourier coefficients 


a 0 


1 fL I 

"IT Jo 


g(z)dz 


( 8 ) 


and 


2 fi-l / «7TZ\ 

"■■r: Jo * Wcos lTrr ■ (9> 


Thus, after time r, at T, with diffusivity D, , 

C 0 /, 2 C 0 . /nw/A /nrrz 

c '< z ’''>=Tr + ,?, — 


Xexp 


m t 


( 10 ) 


Next we assume that during the transition from T x to a 
r 2 , the diffusive redistribution is insignificant. In view of 
he high ratio of thermal to solutal diffusivity of the mate- 
rials under consideration, this is readily achievable in 
practice with rapid temperature ramping. Furthermore, we 
ignore the diffusant redistribution due to the thermal ex- 
pansion or contraction of the sample during the tempera- 
ture change. With typical thermal expansion coefficients 
of liquids, this simplification introduces significant errors 
only for simulation temperature changes in excess of 200 °C. 
Based on these assumptions, the initial condition at T 2 
is taken as 


g2U) = C l( z » r l)- 


(ID 


Thus, accounting for the new sample length L 2 , Eqs. (7) -(9) 
yield for the evolution of the concentration distribution at T 2 
[ 10 ] 


C 2 (z,f ) “ 2 a m 2 cos 


m = I \ ^2 


m TTZ 


exp 


m rr 


D^t 


(12) 


where 



FIG. 3. Definition of sample-collimator-detector geometry for 
the view factor calculation. 


a 02“ 


00 


^+2 
Li n — \ 


An 

L 2 



(13) 


with 



nirlA 

ir) exp 



(14) 


and 


2C 0 /, 

° m 2 tmrLi 


CD 

sin(m7r)+ 2 

n = 1 



sin 

Kir 

SH 

7 TL 2 


n 

m 




II 
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Based on this formulation, we have simulated the evolution 
of the concentration distribution during entire experiment 
scenarios, including the multistep ramps to the first and con- 
secutive measurement tempeatures [10]. 


B. View factor, absorption of radiation, and detector noise 

For realistic estimates of the radiation received at the 
detector, the volume of sample visible through the col- 
limation hole must be known. Hence we have calculated 
the view factor of the sample-collimator-detector geo- 
metry [10]. Assuming perfect absorption by the shield 
material, radiation escapes only through the collimation 
bore of diameter 2 r c and length L c \ see Fig. 3. Under 
this assumption, the number of emitters “seen” by the 
detector is 


J 


N(z d ) = 


r*n f r c rz d + 
J - -nfl Jo Jz d ~i 


* 0 
0 



f I h\ 

./Ml 


r 2 ' 
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/ 1 = 


. — (\/r 2 cos 2 /3+Zo-O 


2 (/,+ r sin /3) 

X 0( \/r 2 cos 2 /3 + Zq _ r <•)• 


(17) 


Zo 


= [L c + 2(/,+ r sin /3)]j cos 2 (18) 

r = (L + / ) and C(z) at the time of measurement is given 
by an equation similar to Eq. (12). As defined in Fig. 3. 
r B,z are the coordinates within the sample, with z at e 
end’ containing the radiation source of initial thickness , . 
The angle /?= - tt/ 2 is taken in the direction of the detector 
on the collimator axis. 

In the actual simulation runs, we have included the attenu- 
ation of the radiation by both the emitter-location-dependent 
self-absorption (see Fig. 2) and the thermal radiation shields 
(see Fig. 6). Furthermore, the noise in the detectors, i.e., me 
uncertainty in the detector signals due to radiation counting 
statistics, was accounted for. This noise is of the order of yfn 
where n is the number of counts received. Thus the statistical 
uncertainty increases with counts, but the relative uncertainty 
decreases as 1 ifn [9], Hence the higher the count rate the 
more accurate the results. This can be accomplished by either 
using a long counting period or high initial activity. 


c. Sample dimensions and activity 

The lengths of the activated section 1, and the overall 
sample L were optimized based on the following consider- 
ations. The above square-root dependence of the detector 
noise on sample activity suggests the use of high activities to 
maximize the resolution in the difference of the two detector 
signals This in turn minimizes the measurement time re- 
quired at a given temperature and thus the errors from con- 
tinuous diffusion during a data acquisition penod. Similarly, 
the initial activity required to ensure a given detector signal 
increases with sample length. However, low activity levels 
are desirable for radiation safety reasons. Simulation results 
(see Sec III E) showed that an initial activity of 5 mil pro- 
vides sufficient resolution. With this initial activity, the op- 
timum sample length is 30 mm [10]. 

To determine the length of the activated section /; we 
recall our assumption that the part of the concentration pro- 
file viewed by a detector is linear; see Eq. (3). us e 
initially activated section must be short enough that the tran- 
sient in which the spreading concentration distribution pos- 
sesses significant curvature is over when the signal level at 
the second detector is high enough to begin measurements. 
The simulations yielded /, = 1 mm as optimum. Thus our i - 
fusion samples consist initially of two parts: the 1-mm- ong 
radioactive section and a 29-mm-long non-activated section. 


tion escaping through the collimation bores. Furthermore as 
can be seen from Fig. 3, due to the conical shape of the 
radiation beam transmitted by the collimation bore, r e s 
larger than the sample diameter will not result in a commen- 
surate increase of the radiation dose received by the detector j 
These arguments are further complicated by the scattering ■ of 
and fluorescence induced by the primary radiation [9j. Our 
samples ('“"In) emit about 95% y and 5% (3 radiation, with 
the average energy in excess of 750 keV. Both cause fluo^- 
rescence in the shield. In addition, y radiation is Compton 
scattered in the shield and, in particular, upon grazing inci- 
dence on the collimator surface. Each scattering process re- 
duces the energy of the radiation. Thus multiple scattering 
and fluorescence can lead to photons that are counted in en- 
ergy windows different from those corresponding to the pri- 
mary photons and f3 particles. ... 

We have experimentally investigated various colhmato 
geometries for their effect on total radiation receive y t e 
detector. Hence we used lead for the collimator optimization 
experiments that were conducted at room temperature. Col - 
luxation holes of 1,1-5, 2, and 3 mm diameter were drilled 
into lead shields of 12.7, 25.5, and 38.3 mm ihickness^ A 
3-mm-diam cylindrical radioactive mdium sample s 
placed consecutively behind the collimation bores. The de 
lector signals were energy discriminated; see Sec. IV G. 
Table I summarizes the results. The signal-to-noise ratio is 
defined as the signal amplitude obtained in * e e nergy wnv 
dow (±0 3 keV about the chosen energy) divided by the su 
of amplitudes measured in all other channels outside the en- 

^From these data, we can draw the following conclusions. 
Since the signal-to-noise ratio (SNR) for the 1-mm coll.ma_ 
tor hole is essentially 1, the smallest usable diameter is 1 5 
mm. As expected, with the 3-mm sample diameter very^u 
additional radiation is received on expanding the colhmato 
diameter to more than 2 mm and, hence, the SNR remams 
essentially unchanged or even decreases. Funher ^^ * e 
signal-to-noise ratio increases with the increase in the thick- 
ness of the radiation shield, from 12.7 mm to 25.5 mn ^ 
a 38 3 mm thickness results in little further improve - 
Copper initially planned as the high-temperature radia- 
tion shidd material, showed intolerably high fluore ^ nce 
levels. For the actual high-temperature expenments, the ^ ra- 
diation shield consisted of gold with a wall thickness of 5J 
mm As calculated from absorption data [9], this wall thick- 
Tss provides the same absorption for 190-keV photons a 
the 25 5 mm of lead used above. This results in an overall 
diameter of the radiation shield/isothermal liner of 41 mm. 
With this collimator geometry and 10-mm squ^ dejetors, 
the maximum distance between the outside edge of the col- 
limator and detector must not exceed 30 mm if all radiation 
from the collimators is to be collected. 


D. Determination of collimator geometry 

Geometric factors as well as radiation properties of the 
shielding material must be considered in optimizing the col- 
limator geometry. Obviously, the radiation shield must be 
thick enough to absorb a large fraction of the radiation. Oth 
erwi.se the detector signal is not representative of the radia- 


E. Simulated multitemperature experiments 

Based on the above considerations, utilizing Eqs. (1)-(18) 
evaluated for the optimized dimensions of our expenmenta 
hardware, simulations were carried out to exp ore the feasi 
bilitv and accuracy of measurements at several temperatures 
Sr same sample. Table II shows to, selMiffus.on ,» 
indium that a temperature range of 700 °C can be covere 
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TABLE l. Signal-to-noise ratios obtained with various collimator geometries. 


Diameter 

(mm) 

Lead thickness 
(mm) 

Signal to noise ratio 
at 24 keV 

Signal to noise ratio 
at 190 keV 

1 

12.7 

1.1 

1.1 

1.5 

12.7 

1.4 

1.25 

2 

12.7 

1.67 

1.5 

3 

12.7 

2.1 

2.0 

1 

25.5 

I 

1 

1.5 

25.5 

1.6 

1.6 

2 

25.5 

2.0 

1.8 

3 

25.5 

2.0 

2.1 

1 

38.3 

1 

1 

1.5 

38.3 

1.9 

1.5 

2 

38.3 

2.2 

1.9 

3 

38.3 

1.7 

2.2 


100 °C increments using only two samples. The second col- 
umn lists the input diffusivities [13] used in the simulations. 
The table illustrates well the importance of taking the detec- 
tor noise into account. Without superposition of noise on the 
signals n { (t) and n 2 (t) y the simulations unrealistically yield 
an error in all eight deduced diffusivity values of less than 
0.2%; see the third and fourth columns. However, as shown 
by columns six and seven, on superposition of a realistic 
±(0)2 noise level, the errors increase to (a still very 
acceptable) 1.5%. 

Note that it is advantageous to begin with the highest 
temperature and, thus, highest diffusivity in a measurement 
series and then to take data at consecutively lower tempera- 
tures. Obviously, at the highest diffusivity the time required 
until the signal at the second detector rises above background 
is minimized; see Eqs. (7)-(9). In addition, the simulations 
unexpectedly revealed another advantage of decreasing tem- 
perature sequences, in which the measurement time at each 
temperature is shorter than at the immediately following one. 
While at all temperatures the spreading of the concentration 
profile depends on the product of the respective D and mea- 
surement time [see Eqs. (7) and (12)], beginning with the 
second temperature there is an additional time-dependent 
spreading term [see Eq. (14)]. As a consequence, longer 
measurement times at the prior temperature, as they would 


be required in a sequence with increasing temperatures, 
cause a larger decrease in the overall concentration gradient. 
The simulations show that the resulting smaller differences 
between the two signals lead to fewer feasible temperatures 
in an experiment. 

The above simulations provided valuable guidance for ex- 
periment design. However, they do not account for some of 
the error sources present in actual experiments. These in- 
clude imperfect shielding and secondary scatter by the radia- 
tion shield and collimation bores. In Sec. IV D we will 
present procedures to at least partly correct for these errors 
using signals from background detectors. There are also sig- 
nificant systematic errors associated with the dimensional 
tolerances of the sample and collimator-detector geometry. 
These error sources are evaluated in Sec. IV E. 

IV. EXPERIMENT 

A. Sample preparation and containment 

The sample consists of two parts: a short radioactive sec- 
tion and a long nonactivated section. The preparation and 
containment of these sections is based on the following con- 
siderations. Uncertainties in the sample geometry are the ma- 
jor error source in these experiments. Furthermore, the 
sample cross section must be uniform to assure one- 


TABLE II. Comparison of simulation model input diffusivities and simulation results for two model runs 
at eight different temperatures. 


Temperature 

(°C) 

Input 

10 5 Z) (cm 2 /sec) 

Result without noise 
10 5 D (cm 2 /sec) % error 

Result with noise 
10 5 £> (cm 2 /sec) % error range 

900 

10.2226 

10.2002 

0.22 

10.2589 

<0.5 

700 

8.260 01 

8.255 02 

0.060 

8.185 63 

<1 

500 

5.978 81 

5.978 31 

0.0084 

5.956 64 

<1 

300 

3.456 23 

3.456 21 

0.000 72 

3.498 95 

<1.5 

800 

9.280 63 

9.25 65 

0.26 

9.331 37 

<0.5 

600 

7.158 44 

7.154 24 

0.059 

7.100 53 

<1 

400 

4.734 12 

4.733 54 

0.012 

4.723 77 

<1 

200 

2.210 22 

2.2102 

0.000 83 

2.210 52 

<1.5 
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FIG. 4. Steps for preparation of the enriched isotope pellet for 
activation. 


dimensional diffusion. Nonuniformities in the cross section 
can arise from capillarity (rounded edges at ends of both 
sections) and voids or bubbles in the sample bulk and at the 
interface between the two sample sections. Equally impor- 
tant is the accuracy in sample length and volume; see Sec. 

IV E. In addition, oxidation, in particular of the contacting 
ends of the sections, can form a diffusion-resistant layer. 
Hence, both the indium and boron nitride (BN) (Carborun- 
dum grade AX05) container parts (BN is not wet by indium) 
are thoroughly high- vacuum baked (10 mm of Hg at 
900 °C) to remove gaseous contaminants. 

The radioactive section is prepared, as shown in Fig. 4, 
using a boron nitride crucible consisting of three parts ma- 
chined with tight tolerances. After insertion of the bottom 
plug into the 3-mm-i.d. cylinder, a bead of indium (40% 
enriched in U3 In) is put in place. The weight of the bead is 
chosen such that the thickness of the final disk will be ap- 
proximately 1 mm. After insertion of the top plug, the whole 
crucible, with a small weight on the top, is inserted into a 
glass tube. This is followed by vacuum heating to 900 c 
and solidification of the sample under high vacuum while it 
is compressed by the weight. After cooling, we backfill with 
argon. Then the section is prepared for irradiation. Utilizing 
a glove box, the glass tube is broken and the sample is trans- 
ferred to a 3-mm-diam glass capillary. To facilitate sample 
removal after irradiation, 5-mm-long rods are placed above 
and below the sample prior to sealing of the capillary. The 
sealed-in section is then sent to the reactor for activation by 
fast neutron irradiation. 

The 29-mm-long nonactivated section of the sample is 
cast in place under high vacuum at 900 °C in the boron ni- 
tride ampoule used for the diffusion measurements. As de- 
picted in Fig. 5, this 7-mm-o.d. ampoule is 60 mm long and 
has a 3-mm precision-reamed bore. To ensure the proper 
sample length we require square comers on the top and bot- 
tom ends of the sample. At the bottom this is achieved by 
press fitting a BN plug into the ampoule bore. At the top we 
apply a slight pressure on the sample with a weak graphite 
spring [Energy Science Laboratories, force 29.4 mN (3 g) at 
900 °C] against a BN plunger that has less than 0.025-mm 
clearance in the bore. To establish conditions that ensure the 
absence of void formation, inactive sections were repeatedly 
melted, solidified, and x-ray inspected. 

Then the activated section is transferred into the ampoule 
in a glove box flushed with an inert gas. After replacing the 



plunger and spring, the BN ampoule is mechanically closed 
with the BN screw lid and inserted into the square-bottom 
silica tube, which,. in turn, is evacuated and sealed. 


B. Diffusion measurement apparatus 

The apparatus, which accommodates the above sealed 
liffusion sample cartridge, consists of the following com- 
wnents (see also Fig. 6). There is a gold radiation shield 
md isothermal liner (90 mm long, 41 mm in diameter)^ 
3old was chosen because of its high melting point and 
hermal conductivity together with the low fluorescence 
md high radiation absorption properties. To minimize dis- 
placement of this gold cell by thermal expansion it is sup- 
ported by a silica pedestal (thermal expansion coefficient 
5.5X10" 7 cm/cm °C). There is a stainless-steel-encased, 
ceramically insulated resistance heater (Fast Heat, 130 mm 
long 120 V 5 A). To avoid radiation absorption, the heater 
has six holes (7 mm in diameter) at the positions of the 
collimators and background detectors. Thermocouples (type 
K) are inserted into the gold liner for temperature uniformity 
control One of these provides input to a temperature con- 
troller and silicon controlled rectifier that stabilizes the tem- 
perature to better than ±0.5 °C. An inner thermal radiation 
shield consists of a stainless-steel sheet. An outer therma 
shield consists of gold-coated boro-silicate glass. There are 
six radiation detectors (not shown): one pair each aligned 
with the oppositely placed, differently spaced radiation bores 
(see Sec. II B) and one pair 90° off the measurement-detector 
plane to record the background radiation. Top and base 
plates contain the measurement structure. Axial heat losses 
from the heated components to these plates are minimized by 
Fiberfrax insulation. 




57 


REAL-TIME D1FFUSIVITY MEASUREMENTS IN . 


1733 


Fiberfrax Top plate 



0 3cm 


FIG. 6. Cross section of diffusion measurement apparatus (with- 
out detectors and electronics). Note the two pairs of collimation 
bores, offset with respect to each other, to cover the temperature 
range 300 °C-900 °C for indium. 

The temperature uniformity across the diffusion sample in 
this apparatus was checked using a dummy diffusion sample 
cartridge with eight built-in, calibrated thermocouples (type 
K y sheath diameter 0.5 mm. Omega Engrg.). At 500 °C, the 
vertical temperature difference A T v across the indium 
sample (see Fig. 5) was 0.4 °C, with the bottom cooler than 
the top. At the same time the maximum horizontal A T hl 
which occurred at about midheight of the sample, was 
0.3 °C. At 900 °C, we obtained a A7 U = 1.4 °C of the same 
orientation and a maximum A7 /t = 0.4 o C at the top of the 
sample. 

An important design consideration for the use of the ap- 
paratus over a wide temperature range is the thermal expan- 
sion of the components, which changes the spacing of the 
collimation bores and their position with respect to the de- 
tectors. The changes in the bore spacing, which over a tem- 
perature range of 700 °C is about 0.4 mm, is taken into ac- 
count when applying Eq. (3). The changing bore locations 
with respect to the fixed detector positions are taken into 
account by using detector crystals large enough that the ema- 
nating radiation cone remains fully detected during the ther- 
mal displacement. 

C. Detectors, discriminators, and calibration 

After examining numerous detector types, we have found 
that CdjZn^Te crystals (eV Products, 1 cm 2 and 3 mm 
thick) provide a quantum efficiency of approximately 1 for 
the 24-keV photons and 0.35 for the 190-keV photons of 
indium. In addition, these detectors exhibit a low noise 
threshold (^10 keV) at temperatures ^50 °C, i.e., no addi- 
tional cooling is required. Two single-channel analyzers 


(SC As) per detector can be set using four potentiometers on 
the circuit boards providing a 0-250 keV energy range. 

Due to differences in the characteristics of individual 
crystals and detector circuits each detector must be indepen- 
dently calibrated. The calibration is accomplished using a 
multiple-channel analyzer (MCA) (Accuspec, Nal, Canberra 
Instruments) and a data acquisition system with a counter tor 
each energy channel. The calibration of each detector SC A 
energy window is performed by matching the SCA output 
count rate and the chosen MCA energy window count rate. 

D. Measurement procedure 

After loading the sealed sample ampoule into the gold cell 
and repacking the thermal insulation, the temperature is 
ramped at 15 °C/min to the highest experiment temperature. 
After the signal at the 5L/6 detectors has sufficiently risen 
above background, data acquisition begins. The measure- 
ment times at this first and consecutive, lower temperatures 
are predetermined based on the simulation results. After dif- 
fusant uniformity has been obtained throughout the sample, 
we return to each experiment temperature and collect base- 
line averages. These data are used to correct for inconsisten- 
cies of the detector crystals and possible abnormalities in the 
collimator-detector geometry. One detector-channel is se- 
lected as being ‘ ‘correct’ * and, using the averages at each 
temperature, all other channels are proportionally adjusted. 
The same is done with the two background detectors. The 
background is subtracted from all channels and then Eq. (3) 
is applied to find the concentration at all times. 

With sources activated to 5-10 mCi we can typically ob- 
tain data at three separate temperatures per run. The simula- 
tions suggested the feasibility of four temperatures. How- 
ever, in the simulation we did not allow for cool down 
periods between temperatures. At lower temperatures these 
periods can be up to an hour for a 200 °C change. 

E. Estimation of uncertainty in apparent diffusivity 

As expressed in Eq. (2), the diffusivity D is obtained from 
the slope of a straight line fit to ln[n i(f) ^(Ol- We have 
shown earlier that such straight line plots can be obtained 
even in diffusion experiments with considerable convective 
contamination [14]. Hence, considering the deviations from 
isothermality of our samples (see Sec. IV B) and thus likely 
resulting convective transport, the following estimate of un- 
certainty in D must be understood with respect to the appar- 
ent rather than the true diffusivity. This uncertainty can be 
expressed in terms of the uncertainties in the slope and the 
sample length as [15] 



where the slope B= - tt 2 DIL 2 and <r 2 B = Na 2 y / A, with N the 
number of data points used in the fit, and A = N('Z i=l t i ) 
-(Z?.,/,) 2 and o- 2 = [l/(/V-2)]S?L I (y ( -A-flf I ) 2 . with 
y = ln(n 2 - «i) and A (the y intercept) = const. 

The uncertainty in the length of the sample can be esti- 
mated as follows. The sample is made of two sections (see 
Sec. Ill C) producing an uncertainty 
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TABLE III. Apparent self-diffusivities obtained with indium at 
the iwo photon energies in four experiments. 


Temperature 

(°C) 

190 keV 
(10 5 D cm : /sec) 

24 keV 

(10 5 D cm 2 /sec) 

800 

9.84±0.59 

10.07 ±0.61 

600 

7.69±0.48 

8.02±0.50 

400 

5.25±0.34 

5.01 ±0.34 

700 

7.72±0.47 


500 

6.20±0.39 


300 

4.60±0.33 


900 

8.95±0.54 

8.61 ±0.52 

800 

9.80±0.59 

10.1 6 ±0.6 1 

700 

8.93 ±0.55 

9.62±0.60 

900 

11. 75 ±0.70 

1 1.63±0.59 

600 

8.04±0.49 

8.09±0.50 

300 

4.54±0.29 

4.78 ±0.3 1 


<SL=<5/ ( +(5/ na , (20) 


where Z NA is the length of nonactivated sample. The uncer- 
tainty in'the two sections is determined from the volume and 
machining of the crucibles such that 



FIG. 8. Plot of 600 °C data of 900 °C-300 °C according to Eq. 
(2), together with diffusivities obtained from least-squares fits 
(straight lines). 

Substituting Eq. (22) into Eq. (21) and the result into Eq. 
(20), the uncertainty in the length of the sample is ± 1.02 mm 
for indium. 

The <j b error is directly related to the counting statistics. 
This error is minimized with high initial activities. With the 
activity levels chosen for these experiments we find that the 
second term in Eq. (19) is up to an order of magnitude larger 
than the first. 



with r the radius of the crucible and V the volume of sample 
material. If we assume that the density p is accurately 
known, then the uncertainty in the volume is 


5V = 



( 22 ) 


The tolerance on the crucibles is ±0.0127 mm (nominal 
sample mner diameter 3.01 mm). Due to the repeated 
vacuum melting, the uncertainty in the masses is 0.001 g. 



FIG. 7. Original signal traces for 190-keV photons during the 
900 °C-300 °C run. 


V, RESULTS AND DISCUSSION 

Several diffusion runs with indium were carried out. 
Each run covered three consecutively lower temperatures. 
The results are summarized in Table III. The error esti- 
mates are based on Eq. (19). The lack of 24-keV data for 
the one run is due to failure of one low-energy channel. Note 
that, within the uncertainty of the data, no difference be- 
tween transport (diffusion) in the bulk and near the wall is 
apparent. 

As-measured time traces of the 190-keV signals obtained 
in the 900 °C-300 °C run (Table III) are shown in Fig. 7. 
These plots reveal several features of the method. The right 
top detector has the highest initial count rate since its posi- 
tion z i — L/6 has been chosen for the highest T= 900 0 C, i.e., 
the longest L\ see Sec. II B. The change in the diffusion 
coefficient at the three temperatures is visible in the slope 
changes of the data. Furthermore, the contraction of the 
sample on temperature decrease is evident from the rise in 
the count rate. One sees also that the final count rates of the 
detectors differ, due to the individual characteristics of the 
detector crystals; see Sec. IV D. A plot of ln(n | / 12 ) vs ti me > 
resulting from similar curves obtained at 600 °C after their 
base-line corrections, is shown in Fig. 8 together with the 
least-squares fits for both energies. 

All data from Table III together with the self-diffusivities 
of In measured by other authors [4,13,16,17,18], are plotted 
in Fig. 9. The wide scatter of the data illustrates well that 
diffusivities obtained in liquids at normal gravity are prone 
to be contaminated by uncontrollable convection. As empha- 
sized for liquid diffusivity measurements by Verhoeven [19], 
any horizontal component of a density gradient results in 
convection without a threshold. We have recently estimated 
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FIG. 9. Experimental data for self-diffusivities in indium: O this work; ▲ Foster and Reynik [17]; A Lodding [13]; • Mathiak et al 
[16]; V and V Carreri, Paoletti, and Vincentini [4] with 1.6- and 0.83-mm-diam capillaries, respectively. Ogloblya, Lozovoi, and Chumakov 
[18] with 0.5-mm-diam capillaries. 


convective contributions to transport in our experiments 
[14], The numerical modeling results for self-diffusion in 
In at 730 °C in 30-mm-long capillaries, 3 mm in diameter, 
show that at normal gravity convective contributions be- 
come significant when an applied horizontal temperature 
non-uniformity A T h exceeds 0.01 K. Often, flow due to 
horizontal density gradients can be reduced by keeping, as 
is the case in our experiments (see Sec. IV B) the top of 
the sample slightly warmer than its bottom. However, the 
above simulations [14] show that this can be rather ineffec- 
tive and under certain conditions may even increase the ap- 
parent diffusivity. 

From the above, it is not surprising that our data show 
evidence of convection in the samples. For instance, the ap- 
parent D at 700 °C deduced from the 900 °C-700 °C run is 
significantly higher than that of the 700 °C-300 °C run. This 
may be due to the longer time the first sample was at high 
temperature before the measurement at 700 °C was made. 
Thus there was more time for convective contamination of 
the diffusive transport. Note, however, that the uncertainties 
for the two measurements at 900 °C do not overlap. We as- 
sign this to possible inconsistencies in the sample prepara- 
tion and convective contributions. The low value may be the 
result of partial blockage of the diffusion path by either a 
void or an oxidation layer between the active and inactive 
sections of the sample. 


IV, SUMMARY AND CONCLUSIONS 

We have developed and systematically optimized a tech- 
nique for the in situ measurement of diffusivities in liquids at 
several temperatures with one sample. This approach circum- 
vents the solidification of the diffusion sample required by 
other methods, which can result in significant errors. In ex- 
periments with ll4m In/In, apparent self-diffusivities were ob- 
tained between 300 °C and 900 °C with an uncertainty of 
±5%. By utilizing the different self-absorption characteris- 
tics of the 24- and 190-keV photons of ll4m In, transport in 
the bulk of the sample and near the container wall was in- 
vestigated independently. No differences were found. There 
is clear evidence for convective contamination of our data. 
To obtain increased accuracy in the diffusivity values, we 
plan to perform these experiments in a low-gravity environ- 
ment. 
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Abstract. In liquid diffusion experiments, the diffusion path can be partially obstructed by 
bubbles/voids, or oxide layers between the segments of a diffusion couple. In addition, there 
have been claims of macroscopic manifestations of 'wall effects’, i.e. a perceived dependence 
of the diffusivity D on distance from the container wall over macroscopic dimensions. We 
have numerically simulated the evolution of the concentration field in a 2D diffusion sample 
in response to such deviations from ID transport. We found that, owing to the smoothing 
of the concentration distribution by radial diffusion, significant deviations from the ideal ID 
concentration field occur only in the axial direction. On evaluation of the concentration field for 
an apparent diffusivity D a , either by the semi-infinite capillary methodology or by the method 
of Codastefano el al (1977 Rev. Sci. Instrum. 48 1650), we found that local obstructions 
of the transport path have to be in excess of half the transport path’s cross-section to result 
in D„ ^ 0.95 D Furthermore, we conclude that the radial concentration gradients found in 
solidified diffusion samples are no evidence for a ‘wall effect’ but are likely to indicate convective 
transport contamination. 


1. Introduction 

Bubbles or voids in liquid diffusion samples [1,2] can potentially result in erroneous 
diffusivity determinations. It may also be difficult to prevent the formation of oxide layers, 
which represent a barrier in the diffusion path, between the segments of a metallic diffusion 
couple Other perturbing phenomena observed in diffusion measurements have been 
genetically assigned to ‘wall effects’ [3-5]. By analogy with surface and grain boundary 
diffusion, it is speculated that the diffusivity D varies over macroscopic dimensions from 
the container wall. However, as pointed out by Nachtrieb [5], irrespective of the possible 
underlying mechanism, it is, of course, not credible that a true wall effect could extend much 
further than a few atomic diameters from the capillary wall, and a reasonable explanation 
in terms of an experimental artifact must be sought. 

In the following we will examine to what extent such deviations from ID transport, 
whether real or conceived, influence the evolution of the concentration distribution in a 
diffusion sample. In addition we will evaluate the resulting non-ideal distributions in terms 
of apparent diffusivities D a utilizing both the semi-infinite sample methodology [6] and the 
method of Codastefano el al [7]. 

In the semi-infinite technique [6], the diffusion sample is assumed to be infinite in one 
direction, i.e. with a diffusion path of length L, the initial thickness h of the diffusant source 
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(W,,h uniform concentration C„) must fulfil Ihc corn), non /, « L. This condttion results in 
the concentration profile at time / 

Co 7 ' 2 


C(x,f) = 


: exp 


(£)• 


( 1 ) 


_ nf in C aeain'st x 2 produces a straight line with the slope given by - l/^Dl). 

Thus* a plo 8 . monitors the solute concentrations against time at 

In the Codastefano technique [7], ^* /#\ resDcctively, 

the ositions x = L/6 and 5L/6. The resulting ume traces C,(0 and C a (l). respect, y, 

are related to the diffusivity through a straight line fit in the form 


ln[C,(f) — C 2 (01 = constant 




( 2 ) 



Figure 1. Geometry and modifications of simulation 
in liquid, (b) barrier at the source-solvent interface, 
indicate regions of initial diffusant concentration Co- 


model. Diffusion capillary 
and (c) y-dependent Ds. 


with (a) voids 
Shaded areas 


2. Simulation model 

The basic geometry of the 2D simulation model [8], we 

simulations are depicted in figure 1. In accor an 3 mm an(J initia i diffusant layer 

chose a total sample length ^ f condition s resulting in deviations from ID 

thickness h = 1 mm. Three diffe ^ ^ e barTier ^ one waU> t0 simulate a 

diffusive transport were assumd^ ( ) A £ these parti cular positions straddle the 

void, was centred at x = 3 or x - / mtm t planar bamer 

Codastefano measurement coordinate L/6 at x = 5 mm (figure 1(a);. 
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Figure 2. Normalized concentration profiles at three diffusion times. Full curves, with void at 
x = 3 mm; broken curves, unobstructed diffusion. 


in the source-solvent interface was extended from one wall to either y = 0.5 or 0.75tx> 
(figure 1(b)). (c) A y-dependent diffusivity was assumed, with either D = 0.5 of the input 
diffusivity D 0 at the wall and then a linear rise to D 0 at y = 100 /xm, or D = 10~‘ D 0 at the 
wall rising to D 0 at either 100 or 300 /xm; see the schematic D-profiles 1-3 in figure 1(c). 
The value of 100 /xm represents the range at which a change attributed to the ‘wall effect* 
was determined in [6]. The distance y — 300 /xm corresponds to the sample depth from 
which 24 keV photons are received in our self-diffusivity measurements with indium [8]. 

Since isothermal self-diffusion was assumed, the fluid density is uniform throughout the 
diffusion process. Hence, solutions of the 2D diffusion equation 



where C is the solute mass concentration and i = 1, 2, fully describe transport in this case. 
We used Do = 10.22 x 10 -5 cm 2 s" 1 in all our numerical simulations, which corresponds 
to the self-diffusivity of indium at 900 °C [9]. 

The numerical simulations were performed with Adaptive Research’s CFD 2000 finite 
volume code, version 3.03. In the void and barrier cases a non-uniform grid of 62 x 30 
and 61 x 40, respectively, was used. For the wall investigations only half the cell 
(0 ^ y < 1.5 mm) was modelled, utilizing a 100 x 61 grid. Doubling of the grid number 
in the x- and y-direction in both the y = 3 mm void and 75% barrier cases (the cases with 
largest concentration gradients) resulted in insignificant changes in the concentration field. 
All cases were run for 10 000 s to obtain sufficient diffusant concentrations at the 5L/6 
position. For comparison with pure ID diffusion an unobstructed, uniform-diffusivity run 
was performed with each grid. 

Diffusivities are typically deduced from the total amount of diffusant contained in cross- 
sectional slices of the sample. Correspondingly, we base our evaluations of the apparent 
diffusivity D a in the various cases on axial profiles of the concentration summed across y, 
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C(.v, t) = C(x, y, t). 
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(4) 
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3. Results and discussion 

Figure 2 shows C(x)-profiles obtained in the case with a void at * = 3 mm a , ,h 

in thTrD‘rffLoMv C oTfreeTL7 loll sets ^ 
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Table I. Apparent diffusivities (JO -5 cm 2 

s ) obtained in void configurations. 



ID transport 

Void @ 3 mm 

Void @ 7 mm 

Semi-infinite method 
Codastefano method 

10.44 ± 0.174 
10.2000 ± 0.0001 

10.18 ±0.72 
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9.9800 ± 0.0005 


■^rSSJS 2£T££ tZZZZZZ‘1 ,D r' 

ID transport run, evaluated with the sem, i„fi„ 7 ! fitS ' Inter estmgly, t 

diffusivity compared wnh ,he inpj * T]o 22 x !JTS ” 

rr No ,e , h a, ~ « , 

sszszzzz ; e snr ° f *- — ■— «- * 

h “ cff “ K - ,i 

is affeced e qM „y regards rf 

increase/decrease in D a with the void upstream/downstrea^ of 2 , / fi S * 

location. This, plus the drastically, lower standard h 7 7 ,77 /6 measuremer 

indicated in tables 1-3 results from rh P . tandard devia,lon of lhe second methodolog 
procedure „U1M ‘he W pro e . 7 i " P, “ ““ ** fol fi “» 

<■» » on a finite resolo, L grid wift 62 pS ' ’-Zc^JT"^ nL““ V ° id ’ “ 
hand, benefits from (a) the smooth variation 7f r/ 7? ^odastefano method - the othe 
and 5Z./6, i.e. jc = 5 and 95 in fi o\ °" °/ C ^ at the two sampling locations ( L/t 

Figure 4 shows the axial profit CM/C dS? PO ' ntS ^ fit ‘ 

at three diffusion times together with thi l ° W ‘^ the 75% barr,er (figure Hb)] 

from the unobstructed ID diffusion case rwT" Ttf ^77 reSUlt ' ng on the same gnd 

a,irus,on case - One sees that this barrier causes a pronounced 


Semi-in 

Codaste 


retenti 
that c\ 
above 
more * 
Th 
the pa 
uncert 



On the insensitivity of liquid diffusivity 
0 


measurements 


7117 


(4) 

* 

0.5 1 

profile 
o et al 


1.0 

ing on 


1.5 



2.0 - 


tf 



£ 

2.5, 


■ three 

X 

times 

3.0 

iect to 


action, 

3.5 

by the 

4.0 ; 

yre 3. 


raight 

4.5 ^ 

acting 



5.0 - 



- 1 i L_ 


0 0.5 1.0 1.5 2.0 2.5 3.0 

y [mm] 


Figure 3. Normalized isoconcentration lines Ctx vWr , u 

t = 1000 s with void at * = 3 mm. > /C " ,n ,he u PP er P 3 * of the capillary at 


suits, 
the 
fated 
izing 
« L 
nt of 
ental 

i the 
uque 
s an 
ffient 

logy 

tting 

and 

*ther 

i/6 

fit. 

(b)) 

grid 

iced 


Table 2. Apparent diffusivit.es (ICT 3 cm 2 s -') obtained in 


I in barrier configurations. 

ID transport 50% barrier 75% bamT 

Semi-infinite method 10.45 ± 0.19 (0 05 ± 0 42 

Codastefano method 10.2000 d: 0.0002 
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Figure 4. Normalized concentration profiles at three diffusion times. Full curves, with 75% 
barrier at x = l mm; broken curves, unobstructed diffusion. 


Table 4. Maximum normalized radial concentration variation at x = Lj 6. 


Case 

Radial AC/C (%) 

Time of occurrence (s) 

D(>) case 1 

0.014 

250 

D{y) case 2 

0.028 

230 

D(y ) case 3 

0.082 

230 


The insensitivity of the concentration distribution to postulated ‘wall effects’ is further 
illustrated in figure 5, in which we have plotted the concentration difference (normalized 
to the samples maximum C) between case 3 and the ID transport result at 300 s. It is 
evident that solute redistribution by radial diffusion prevents the formation of significant 
concentration gradients across the width of the sample. This is also supported by the low 
values of the maximum radial concentration variation at position x - L / 6 listed in table 4. 

Note that these maxima occur early in the diffusion runs. Thus they are unlikely to be 
detectable in diffusivity measurements. As a consequence, experimental findings of radial 
concentration gradients in cylindrical diffusion samples indicate convective contamination 
(for detailed examples see [12]) rather than inferred dependences of diffusivity on the 
distance from the container wall. 


4. Conclusions 

The above 2D simulations show that the influence of the model voids, barriers and ‘wall 
effects on concentration distributions in diffusion samples and the subsequently deduced 
diffusivities is rather marginal. Note that in real, 3D systems, owing to the lower cross- 
sectional fraction taken on by such obstacles of the same radial dimensions, the deviations 
from unobstructed diffusive transport will be even less significant. Hence, we conclude that 
the contamination of diffusivity measurements owing to such geometrical parameters will in 
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Figure 5. Difference in concentration between the ‘wall effect’ case 3 and uniform diffusion, 
A C(x, y), normalized to the maximum concentration, obtained in the lop part of the sample at 
/ = 300 s. Numbers in 10~ 2 %. 


most cases be negligible compared with the convective contamination arising from sample 
non-isothermality [12]. 
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J . I. D. Alexander , J.-F. Ramus , and F. Rosenberger 

Numerical Simulations of the 
Convective Contamination of Diffusivity 
Measurements in Liquids 


Ideally y self ’diffusivity measurements should be conducted 
under isothermal conditions. This is difficult to achieve in 
practice , particularly at high temperatures. Hence , we have 
assessed the effect of buoyancy -driven convection in liquid - 
filled vertical diffusion capillaries with slight azimuthal ther- 
mal asymmetries on their walls and ends. A tracer of same 
density as the liquid is initially located at one end of the 
capillary. Calculations of the 3D -time -dependent transport of 
the tracery based on properties of liquid indium at 1,000 Ky 
yielded a variety of results with practical implications. We 
found that with horizontal temperature nonuniformities 
across the sample as low as 1 and 0.1 K, respectively convec- 
tive transport in capillaries of 1 and 3 mm diameter can 
exceed diffusive transport. Superimposed vertically “ stabiliz- 
ing ” temperature gradient did not sufficiently reduce these 
convective transport contributions. Hence , the widely used 
approach of keeping the top of the capillary slightly warmer 
than the bottom appears to be an ineffective means of “ pre- 
venting convection Most strikingly , we found that diffusive- 
convective transport can result in axial concentration profiles 
that are analytically indistinguishable from pure diffusion 
profiles. This casts considerable doubt on various graphical 
criteria previously used to demonstrate the absence of signifi- 
cant convective contamination in diffusivity data. Further- 
more , these results reinforce the necessity of conducting 
diffusion measurements in liquids under microgravity condi- 
tions where convective contamination can be minimized or 
avoided entirely. 


1 Introduction 

Measurements of diffusion coefficients in liquids by macro- 
scopic techniques are prone to be contaminated by convec- 
tion. For example, in a system with a diffusivity of 
10 -5 cm 2 /s and a diffusion distance of 1 cm, the characteris- 
tic diffusion velocity is of the order of I0~ 5 cm/s. Thus, if 
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convective transport rates are not to exceed 1 % of the 
diffusive transport rates, convective flow velocities must not 
exceed I0 _7 cm/s. In liquid self-diffusion experiments, 
buoyancy arises only due to density gradients introduced by 
temperature nonuniformities. Around room temperature, 
sufficient sample isothermality can be maintained to avoid 
significant convective flows. However, at the elevated tem- 
peratures characteristic of liquid metals and semiconduc- 
tors, it is difficult to realize isothermal conditions to an 
extent that convection remains insignificant. This was rec- 
ognized, in particular, by Verhoevetty who emphasized that 
any horizontal component of a density gradient results in 
spontaneous convection with no threshold [1]. Since then, 
in an attempt to reduce the flow from horizontal density 
gradients by applying a vertical, stabilizing one, it has 
become widespread practice to keep the top of diffusion 
capillaries slightly warmer than their bottom end. 

Convection can be effectively reduced by reducing the 
characteristic distance for momentum transport in the (dif- 
fusing) liquid. However, a systematic experimental investi- 
gation has revealed convective contributions in vertical 
capillaries with diameters as small as 0.7 mm [2]. The same 
work has shown that concentration plots versus the square 
of the distance along the diffusion capillary can result in 
straight lines (often used as criterion for negligible convec- 
tion in semi-infinite-sample techniques) even in the presence 
of pronounced convection. Many workers have attempted 
to detect convective transport contributions by varying the 
diameter of their diffusion capillary. A plateau in the result- 
ing diffusion coefficient versus capillary diameter is then 
interpreted as precluding convective masking; see, e.g., [3]. 
Of course, this conclusion may be erroneous since such a 
plateau only indicates that the convective effect is diameter- 
independent. Other workers have proposed the use of mag- 
netic fields to suppress convection during diffusion studies 
in conducting liquids. Unfortunately, at the very low veloc- 
ities in question, the magnetic field strength required to 
achieve “pure” diffusion is so high that, as other work has 
shown, the diffusion process itself is significantly modified 
through the Lorentz force acting on the ions and electrons 
in the liquid [4, 5]. 

As a consequence of these difficulties, experimental data 
for self-diffusivities of liquid metals and semimetals ob- 
tained with different techniques and/or by different authors 
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Fig. I. Mode! diffusion capillary. 
Tracer initially located in shaded layer 
of thickness S 


typically differ by several 10%, and sometimes by more 
than 100% [6, 7). This is not only an obstacle to realistic 
process modeling, but to our theoretical understanding of 
diffusion in liquids as well. Although various theories exist 
for the diffusion in liquid elements and its tempera- 
ture dependence, the uncertainties associated with most 
diffusivity data are too large to allow for an unambiguous 
evaluation of these theoretical models. Currently, the differ- 
ences between several theoretical predictions are often less 
than those between different sets of data for the same 
system. 

To provide quantitative guidance for experimental work, 
we have numerically modeled diffusive-convective transport 
in vertical cylinders with dimensions typical for diffusion 
experiment geometries. To account for thermal convective 
transport contributions we imposed slight azimuthal tem- 
perature asymmetries on the walls and ends of these capil- 
laries. A diffusant, representing the isotopic tracer used in 
the self-diffusion experiments, is initially located at one end 
of the capillary. The 3D-time-dependent transport of the 
diffusant was then simulated for physical properties corre- 
sponding to liquid indium. In the following, we present the 
simulation model, the numerical approach, and the simula- 
tion results together with their practical consequences. 

2 Simulation Model 

The liquid was assumed to be a Boussinesq fluid contained 
in a closed vertical circular cylinder of length L — 30 mm 
and diameter 2 R (see fig. 1) of 1 or 3 mm. Gravity acts 
along the cylinder axis. The temperature on the vertical 
cylinder wall has the form T(R , 0, z) — (T, + T 2 )/2 + (T ] — 
r 2 )(cos (0))/2 + AT.zjL. Here T X -T 2 = AT ( , is the maxi- 


mum horizontal temperature difference across the cylinder 
at any height z, and AT. is the temperature difference 
between the colder bottom and warmer top faces. The 
temperatures of these end faces are T(t\ 0, 0) = (T, + 7\)/2 
+ (T, - 7 2 )r(cos (0))/(2R) and T(r, 0, L) = 7\r, 0, 0) + 
AT ., respectively. The imposed temperature gradients are 
time-independent. From this, and their assumed small mag- 
nitude, we expect the resulting flow to be steady. Hence, for 
given values of AT 0 , AT., and R we can compute a steady 
solution to the following mass, momentum, and heat trans- 
port equations: 

divtf + 0, 0) 

u • grad u = — V/? T v Au — fi(T— T 2 )gk, (2) 

Q 

u * grad T = k AT. (3) 

Here u, T, <? 0l /?, v, and k are the liquid’s velocity, tempera- 
ture, density, thermal expansion coefficient, kinematic vis- 
cosity, and thermal diffusivity, respectively. Note that, since 
we consider self-diffusion experiments for which the tracer 
density is equal to the density of the host liquid, there are 
no solutal contributions to buoyancy. The no-slip condition 
for velocity is applied at all walls. Having found the steady 
velocity and temperature distributions we then use u to 
compute the evolution of the concentration field c(x\ () 
from 

— + u • grad c — D Ac, (4) 

o/ 

where D is the self-diffusivity. As indicated in fig. 1, the 
diffusant (tracer) is initially located in a thin layer of 
thickness S at one end of the capillary. Hence, the initial 
concentration distribution is 

c(z, 0) = c o , 0 <z<6 and c(z, 0) = 0, S<z<L , (5) 

where S/L « 1 . No flux of the tracer is permitted through 
the walls and ends of the capillary. 

The simulations were carried out using the physical 
properties of molten indium evaluated at T= 1,000 K 
(table 1) together with an assumed dififusjvity of 1-48 ■ 
10“ 5 cm 2 /s. Since it is desirable to measure self-diffusion 
coefficients over a wide temperature range, we have evalu- 
ated the relevant properties also at lower temperatures. One 
sees from table 1 that, in particular, the viscosity varies 
considerably with temperature. For an extrapolation of our 
results to lower temperatures, the reader should multiply 
the r= 1,000 K value of AT 0 by the ratio given in the far 
right column. For example, at 850 K, the convective contri- 



Table 1. Physical properties of indium for different temperatures above the melting point (429 K) calculated from data in [8 9]. The 
simulations were carried out for properties corresponding to T — l ,000 K. To obtain the horizontal temperature difference , A T 0 , which, at the 
lower temperatures would give rise to the same flow as the 1,000 K case , multiply the AT 0 for fOOOK by the ratio given in the last column 


T 

[K] 

e« 

[gem*-’] 

V 

[ 1 0 “ 3 * cm 2 • s _ '] 

K 

[cm 2 ‘ s~‘] 

o 

i 

5 * 

JT„(T) 

AT,, (1,000 K) 

. 1,000 

6.64 

0.78 

0.11 

1.02 

1.0 

850 

6.74 

1.27 

0.11 

1.01 

1 .65 

600 

6.91 

2.09 

0.11 

0.98 

2.79 

429 

7.03 

2.65 

0.11 

0.97 

3.60 
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Fig. 2. Defining sketch for data evaluation; for details see text 


where 6 « L and r(z, t) is bounded as oo. The solution 
has the form [14] 


c(r, /) = 


£o 

V nDt 


exp 



( 10 ) 


and D is obtained graphically from the slope m of 
In (c(z, 0) versus z 2 at t = /*, where m is given by 


flln (c(r, /*)) _ -1 
dz 2 4 Dt*' 


(ID 


In the second method we calculated three diffusivities from 
the concentration profiles along the three vertical axes 
represented by heavy lines in fig. 2. One value was obtained 
from the axial concentration profile at r = 0 (D = D„ y ) and 
the other two from axial profiles taken at r = 0.15R and 
0 = 0 and n (D = D ax0 , D = D axn ) respectively. 


butions observed at T= 1,000 K will occur for a AT 0 value 
which is 1.65 times the AT 0 at 1,000 K. 

The governing equations and boundary conditions for 
the time-dependent transport of heat, mass, and momen- 
tum in the system were solved using finite differences on an 
equally spaced and a variably spaced mesh using the flow 
analysis software CFD2000. With the latter method, spatial 
discretization is achieved using a finite volume technique 
[10]. The PISO algorithm [11] is used to couple the dis- 
cretized heat and momentum transport equations which 
were solved using an ADI method [12]. Owing to the large 
length to diameter ratio of the capillaries (20 and 60) we 
found it necessary to employ 24 x 12 x 140 points in^the 
azimuthal, radial, and axial directions, respectively. 


3 Data Analysis and Presentation 


Time histories of the tracer concentration fields were ob- 
tained for different values of R, AT e , and AT : . To assess 
the consequences of the different temperature nonuniformi- 
ties for diffusion coefficient measurements, we treated each 
of the calculated tracer time histories as we would experi- 
mental data. Two approaches were used to deduce diffusiv- 
ity values from the calculated concentration distributions. 
In the first approach we employed the method of Codas te- 
fano et al. [13] in which time records of the concentration 
c { (t) and c 2 (f), respectively, at the axial sample locations 
z = L / 6 and 5L/6 are used to calculate the diffusivity D 
from the relation 


In [c,(0 - c 2 (r)] = const. ” (n/L) 2 Dt. (6) 

As indicated in fig. 2, the concentration values were taken 
at the respective positions z as averages over either the 
entire cross-section (D = D C ) or over the ring-shaped re- 
gion 0.75 < r < R y 0 < 6 < 2n (D = D w ). 

In the second approach we used the solution to the 
following ID diffusion problem for c(z, t) 


dc _ d 2 c 
Jt~ D d?' 


8c(0,/) _ 

dz 


■■ 0 , / ^ 0 , 


c(z, 0) = c 0 , 0 < z < (5, 


and c(z, 0) = 0, 


(7) 

( 8 ) 

3<z<L , (9) 


4 Results and Discussion 

The results of our simulations are summarized in table 2. 
One sees that, irrespective of the application of a “stabiliz- 
ing” AT with the 1 mm diameter capillary the horizontal 
temperature difference AT 0 must not exceed 0.1 K for the 
apparent diffusivity D to be within a few percent of the 
actual (input) diffusivity D 0 . However, at AT 0 = 1 K, which 
we consider as a readily occurring temperature nonunifor- 
mity for experiments at 1,000 K, D c exceeds D 0 by over 
100%. The magnitude of this convective contamination, 
though somewhat reduced by the A T. of 2 K, fully demon- 
strates the difficulty of accurately measuring diffusivities at 
elevated temperatures in the presence of gravity. These 
results are further illustrated in fig. 3a with time records for 
c, and c 2 used in the Codastefano method, eq. (6). The 
convective transport contributions become particularly ap- 
parent when comparing the curves obtained with gravity 
(heavy solid lines: AT : = 0, dashed lines: zJ7"_ = 2K) with 
the thin solid line resulting under purely diffusive condi- 
tions. On the other hand, the plots of In (c, - c 2 ), obtained 
under these strongly contaminated conditions, are as 
straight as one would expect for pure diffusion. Thus, the 
straight line character of these plots cannot be used as a 


Table 2. Apparent diffusivities in units of [10~ 5 cm 2 s *] obtained 
from the simulations with input diffusivity D 0 = 1.48 ■ 10 5 cm 2 
s -i' Values with * obtained at t = 300 min , with ** at 356 min 


d7>0K 

2 R 

= 1 mm 


2 R 

= 3 mm 


AT 0 [ K] 

1.0 

0.1 

0.01 

0.1 

0.05 

0.01 

D r 

3.69 

1.50 

1.48 

7.92 

3.08 

1.54 

D w 

3.69 

1.50 

1.48 

7.92 

3.08 

1.54 

D ax 

2.25** 

1.50* 

1.51** 

7.87** 

3.26* 

1.60* 

D n 

2.25** 

1.49* 

1.49** 

6.18** 

2.75* 

1.51* 

D ax , 

2.25** 

1.49* 

1.49** 

8.43** 

3.45* 

1.57* 

II 

to 

* 

AT 0 [ K) 

1.0 

0.1 


0.1 

0.077 

0.01 

D r 

2.63 

1.51* 

- 

8.93 

5.73 

1.55 

D w 

2.63 

1.51* 

- 

8.93 

5.73 

1.57 

d' x 

2.25** 

1.48* 

- 

8.76** 

7.02* 

1.60* 

D WJ 

2.25** 

1.47* 

- 

7.19** 

5.96* 

1.51* 

A,,. 

2.25** 

1.47* 

- 

8.76** 

7.64* 

1.57* 
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Time [min] Time [min] 


Fig. 3 (left), (a) Normalized average con- 
centration c, (upper curves) and c 2 (lower 
curves) vs time and (b) In (c t — c 2 ) vs time , 
for 2R = / mm and AT ( ,= I K. The pure 
diffusion case is shown for reference in (a). 
The solid and dashed lines correspond to 
A T. = 0 and 2 K , respectively 

Fig. 5. (right), (a) Normalized average con- 
centration c, (upper curves) and c 2 (lower 
curves) vs time and (b) In (c, — c 2 ) vs time , 
for 2R = 3mm and AT 0 = 0. 1 K. The pure 
diffusion case is shown for reference in (a). 
The solid arid 'dashed lines correspond to 
AT. = 0 and 2 K, respectively. 


criterion for the absence of significant convective contribu- 
tions. More insight on the convective mixing in this case 
can be obtained from the isoconcentration lines in the plane 
0 = 0, 7 c (fig. 4). The deformation of these lines, which 
under purely diffusive conditions would be strictly straight 
and horizontal reflects a convection role with the liquid 
rising at the right (“hotter”) wall and sinking at the left 
wall. Characteristic convection velocities associated with 
these rolls are of the order of 8 • 10" 4 cm/s. A comparison 
of these concentration contours for A T s = 0 (solid line) and 



Fig. 4. Normalized concentration contours for 2R ~ l mm and 
A T„ — IK at (a) t - 178 and (b) 356 min. Contour intervals: (a) 
5- I0- 4 ; (h) 4- I0- 4 


A 71 = 2 K (dashed lines) again illustrates the limited effect 
of the common practice to keep the top of the diffusion 
capillary somewhat warmer in order to suppress convec- 
tion. 

The results for 3 mm diameter capillaries are even more 
sobering. Table 2 reveals that AT ff *s of a few 10“ 2 K suffice 
to induce intolerable convection contamination. For the 
case of AT 0 = OA K, the strong convective contributions 
become apparent by comparing the plots of c x and c 2 in fig. 
5a with those in fig. 3a. One sees that in the wider capillary 
the- concentration at the “Codastefano locations” (see 
above) increases much more rapidly than in the narrower 
capillary. This is due to an increase of the convective 
velocities to 7 • 10 -3 cm/s. This increase is also reflected in 
the stronger deformation of the isoconcentration lines, 
which are presented in fig. 6 for this case. In spite of this 
heavy convective contamination, the linearity of the 
In (c, - c 2 ) versus time curves at large times (fig. 5b) might 
again lead one to erroneously conclude that convection is 
insignificant! Equally interesting is the result that AT. = 
2 K, rather than reducing the flow induced by AT ff , slightly 
enhances the convective contribution. We have observed a 
similar behavior associated with transport during Bridg- 
man -Stockbarger solidification under microgravity condi- 
tions [15]. 

The above discussion was restricted to D -values ob- 
tained by the Codastefano method, eq. (6). The apparent 
diffusivities obtained by the second method, eq. (10), reflect 
a similar behavior. Despite the obvious convective contam- 
inations, plots of the logarithm of the concentration versus 
the square of the distance along the diffusion capillary 
(which are employed in this method for D -evaluations, see 
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Fig. 6. Normalized concentration contours 
for 2R — 3 nun and A T 0 — 0.1 K at (a) f - 
178 and (b) 356 min. Contour intervals: (a) 
10~ (b) 4 • 10~ 4 


Fig. 7 . Plots of In c vs z~ along the three 
axis defined in fig. 2. Same sample as in figs. 
3 and 4: 2R = I mm and AT,, - 1 A, A T. = 
0, at t = 356 min. Dashed straight line: guide 
for the eye to emphasize absence of curva- 
ture in profiles over large parts of the sample 



Fig. 8. Plots as fig. 7 samples as in figs. 5 
and 6: 2R = 3 mm and AT,, — 0. 1 K, A T. = 
0, at t- 356 min 


eq. (11)) result in rather straight lines over major portions 
of the sample. This is depicted in figs. 7 and 8, respectively, 
for the AT : = 0 cases of figs. 3 and 5. Typical liquid 
diffusivity studies involve the solidification of the sample 
before the chemical analysis for c(z) determination. This 
solidification often results in high porosity or deformation 
of the last part of the sample to freeze. On discarding this 
part of the sample, an investigator would again reach the 
erroneous conclusion that convective contamination was 
insignificant. . I 

5 Summary and Conclusion 

We found that temperature nonuniformities of magnitudes 
that, in particular at high temperatures, are practically 
unavoidable during diffusivity measurements in liquids, 
cause strong convective contamination of the diffusivity 
data. Furthermore, the convective contribution to transport 
is only insignificantly reduced by a “stabilizing’ AT z . 
Hence, the widely used approach of keeping the top of 
diffusion samples slightly warmer than the bottom appears 
to be a dubious means of “preventing convection”. Most 
strikingly, we found that diffusive-convective transport can 
result in temporal local concentration changes and axial 
concentration profiles that are analytically indistinguishable 
from pure diffusion profiles. This finding casts considerable 
doubt on various geometrical criteria previously used to 
demonstrate the absence of significant convective contami- 
nation in diffusivity data plots. 

From this study we conclude that reduced gravity condi- 
tions are essential for definitive diffusivity measurements in 
liquid phases of macroscopic dimensions. 
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